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1. Introduction 


The objective of this report is to establish an analytical six degrees of freedom 
(6 DOF) aerodynamic model of a high angle-of-attack (alpha) combat airplane that can be 
utilized in optimization and control analysis/synthesis studies. Emphasis is placed on 
deriving such a model with validity in the altitude-Mach flight envelope centered at an 
altitude h = 15,000 feet and a Mach number M = 0.6. Some effort is made to extend the 
validity from 0.3 to 0.9 Mach. An engine model is not included. The analytical models of 
aerodynamic derivatives are derived as nonlinear functions of alpha with all other states and 
control variables fixed. Consequently, interpolation is required between the parameterized 
nonlinear functions. 

In this report, a six degree of freedom (6 DOF) sub-sonic analytical aerodynamic 
model is derived from a high angle of attack research vehicle wind-tunnel model. The 
wind-tunnel model was provided by NASA LaRC, [1], which is based on that contained in 
[2-3]. The derivation uses only the aerodynamic coefficient data of the wind-tunnel model 
which corresponds to an altitude h = 15,000 feet and Mach numbers ranging from 0.3 to 
0.9. In order to avoid additional complexity, certain effects are not considered: The effects 
of leading edge flap, trailing edge flap, speed brake, landing gear, etc. The aerodynamic 
coefficients are considered to be functions of the following control variables as well as 
angle of attack, sideslip, Mach number, altitude, roll, pitch and yaw rates: Aileron 
deflection (8a), Rudder deflection (8r) and Stabilator deflection (5h ). In addition, lift and 
pitching moment coefficients have unsteady flow parts due to the time rate of change of 
angle of attack (alpha dot). 

Using body axes, the equations of motion are developed in Chapter 2 in which the 
center of mass (eg) and the aerodynamic center(ac) may be non-colocated. The 
aerodynamic coefficients modeled are drag, lift and side forces and rolling, pitching and 
yawing moments. After presenting the mathematical structure of the aerodynamic 
coefficients which has a dependency on alpha dot, explicit equations of motion are 
developed. 

The derived 6 DOF analytical model is presented in Chapters 3-8. The derivation is 
based on a high alpha research vehicle (HARV) wind-tunnel model described in [1]; it is a 
full, nonlinear 6 DOF, rigid-body dynamic model whose aerodynamic forces and moments 
are calculated from a large wind-tunnel-derived data base using table look-ups with linear 
interpolation. The angle of attack range is -10° to 90°; sideslip angle range is -20° to 
20°; Mach number range is 0.2 to 2.0 and altitude range is 0 to 60,000 feet. Only subsonic 
Mach numbers and the fixed 15,000 ft altitude are considered in fitting analytical models to 
the data. 

In Chapter 3 the drag coefficient is modeled at Mach 0.6 using four nonlinear 
functions of alpha which are parameterized by stabilator deflections 8h = 10.5°, 0°, -5° 
and -24°. The formulae are given in Table 3.3 and comparisons with the wind-tunnel data 
are presented in Figure 3.1. In Chapter 4 the lift coefficient is modeled at Mach numbers 
0.6 and 0.9 and parameterized by stabilator deflections 8h = 10.5° and -24°. The formulae 
are given in Table 4.3 and comparisons with the wind-tunnel are presented in Figures 4.1 
and 4.2. 
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Mach numbers 0.3, 0.6, 0.8 and 0.9 are used in Chapter 5 to parameterize the 
analytical models for the pitching moment coefficient. They are also parameterized by 
stabilator deflections 8h = 10.5°, 5°, 2°, 0°, -5°, -12.5° and -24°. The formulae are given 
in Tables 5.3 a,b,c,d and comparisons with the wind-tunnel data are presented in Figures 
5. 1-5.8. 

The analytical model of the side force coefficient C y is given in Chapter 6. It is 
taken from the wind-tunnel model at an altitude h= 15,000 feet and a Mach number M=0.6. 
The analytical model for C y() is constructed at (3 = 0°, 20°; 8a=+ 25°; and 8r = + 30°. The 

analytical models are functions of a from 0° to 90°; they are defined in Tables 6.1 and 6.2. 
The analytical formulae are presented in Table 6.3. Comparisons of the analytical models 
with the corresponding wind-tunnel model data are shown in Figures 6.1 to 6.3. The roll 
and yaw rate derivatives C yp and C yr are given in Figure 6.3. 

The analytical model of the rolling moment coefficient Q is given in Chapter 7. It 
is taken from the wind-tunnel model at an altitude h= 15,000 feet and Mach numbers M=0.6 
and 0.9. The analytical model is parameterized at P = 0°, 20°; 8a=+ 25°; and 8r = + 30°. 
The analytical models are nonlinear functions of a from 0° to 90°; they are defined in 
Tables 7.1 a,b and 7.2 a,b. The analytical formulae are presented in Tables 7.3 a,b. 
Comparisons of the analytical models with the corresponding wind-tunnel model data are 
shown in Figures 7.1 to 7.4. The roll and yaw rate derivatives Q p and Q r and the sideslip 

derivative C^ are given in Figure 7.4. 

The analytical model of the yawing moment coefficient C n is presented in Chapter 
8. It is tiken from the wind-tunnel model at an altitude h=15,000 feet and Mach numbers 
M=0.6 and 0.9. The analytical model is parameterized at P = 0°, 20° and stabilator 
deflections 8h = 10.5° and -24°; 8a=+ 25°; and 8r = + 30°. The analytical models are 
functions of a from 0° to 90°; they are defined in Tables 8.1 a,b,c,d and 8.2 a,b,c,d. The 
analytical formulae are presented in Tables 8.3 a,b,c,d. Comparisons of the analytical 
models with the corresponding wind-tunnel model data are shown in Figures 8.1 to 8.7. 
The roll and yaw rate derivatives C„ p and C„ r and the sideslip derivative C n p are given in 
Figure 8.7. 

The wind-tunnel model of [1] was flown in NASA's simulator by a pilot to 
generate some basic maneuvers at 0.3, 0.6 and 0.9 Mach numbers such as pitch-ups, 360° 
loaded and unloaded rolls, turn reversals, split S's and level turns. That simulator data was 
used to check the validity of the 6 DOF analytical model. The accelerations 

u, w, q, v, p, r 

are computed for the analytical model using the states and controls from the piloted 
simulated maneuvers. Comparisons with the accelerations from the wind-tunnel data 
model are shown in Chapters 9, 10 and 1 1 for the Mach numbers 0.6, 0.9 and 0.3, 
respectively. 

Listings of the computer code developed under this grant are contained in the 
Appendices A-D. Appendix A contains the code for the equations of motion. Appendices 
B and C contain the code for the analytical model. Appendix D contains the code for the 
comparison. 
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2 . Equations of Motions 


2 . 1 Equations of Motions - Body Axes 

The following kinematical relations for a rigid symmetric aircraft are given with reference 
to body axes at the center of gravity, [4]. The state vector is (u,v,w, p,q,r,Q,§). The 
aerodynamic linear acceleration vector at the center of gravity (eg) is denoted as (X,Y,Z). The 
aerodynamic angular acceleration vector at eg is given by (FP.FQ.FR). The thrust vector T is 
represented in body coordinates as (T x ,T y ,T z ). 

The force equations with respect to body axes: 

T 

u = rv-qw- gsind + X + — n 

m W 


v = pw- ru + geos 6 sirup + Y+ — 

m 


w-qu-pv + gcos 9 cos <f) + Z + — 

m 

The moment equations with respect to body axes: 

P = C M pq + C n qr + C n FR + C'FP 

*1 *x 

q = C„ pr + C s2 (r ! - p’ ) + FQ + 1 (/„ T, - *„T, ) 

*y 


r = C 61 pq + C (2 qr + C 62 FP + C'FR 


(4) 

(5) 


T i t ') + ££(£ T -t T ) 

/ m 7 ] ** y ye X 

Y 7 


( 6 ) 
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Euler Equations: 


6 = qcostp-rsin(f) 

0 = p + qtanO sirup + r tand cos <p 


(7) 

( 8 ) 


where he quantities X, Y, Z, FP, FQ and FR depend on the aerodynamic coefficients C D , 

C ,C. , 2,,C m ,C„ as follows: 
y L* c ’ m n 


D = qSC D (Drag) 

(9) 

L = qSC L (Lift) 

(10) 

X = [-Dcos(a) + Lsin(a)]/ m 

(ID 

Y = qSC / m 



y 

(12) 

Z - [-Dsin(a) - Lcos(a)]/ m 

(13) 

FP = 

qSbC , +m{l Z-l Y )]// 



l y y z ) \ x 

(14) 

FQ = 

qScC +mU X-i z)l / / 



m \ z x )\ y 

(15) 

FR = 

qSbC +m(l Y-t X^ // 



n v x y A 2 

(16) 
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and where the constants in the moment equations ( p,qand r) are functions of the moment 
of inertia quantities I X) Iy, I z and I xz : 


2 


C* = l x I//(I X I Z - lx/.) 

07 ) 

C41 = C I XZ (I Z + lx • Iy)/ lx Iz 

( 18 ) 

2 

C 42 =C* (I z (ly-I z )-Ixz)/Ixlz 

( 19 ) 

C43 = C I xz /Ix 

(20) 

C 5 l = (I z -Ix)/Iy 

(21) 

C52 = Ixz/Iy 

(22) 

2 

C 6 1=C* (I x (Ix-Iy) + Ixz)/I x lz 

( 23 ) 

C 6 2 = C* I xz (l y -I z -I x )/I x I/ 

( 24 ) 

c 63 =c* I xz /I z 

( 25 ) 


In the above equations the vector ( £ x , £y. It) denotes the position vector from the 
center of mass (eg) to the aerodynamic center (ac) and the vector (£ x& < fy e - £ Z e) denotes the 
position vector from the center of mass to the engine thrust center. 

The thrust vector T = (T x ,T y ,T z ) has the components: 

T x = T R (h,M,5 T ) cos(1.98°) + T L (h,M, 6 T ) cos(-1.98°) 

T y = T R oi,M,5 T ) sin(1.98°) + T, (h,M,S r ) sin(-1.98°) 

T z = 0 

where 87 is the throttle control variable which varies between 30° (idle) to 131° (full after 
burner). Since only aerodynamic models are considered in this report the engine models 
for computing the thrust, T R (h,M, 8 T), of the right engine and the thrust, T L (h,M, 8 j), of 

the left engine are not provided herein. The thrust angle 1.98° is the position of the engines 
away from the center line of the aircraft. 
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The following constants are typical of those for a high alpha research vehicle 
(HARV), [1], 

Table 2.1 Aircraft Constants for a High Alpha Research Vehicle 

m = 1035.308 slugs 
w = 33310 lbs 
I x = 23000 slugs ft 2 
I y = 151 293 slugs ft 2 
I z = 169045 slugs ft 2 

! xz = -2! >71 
S = 400 ft2 

c = 11.52 ft 
b = 37.42 ft 
= -0.297 

£y = 0. ft 

e 7 = 0.233 ft 
£ xc = -19.37 ft 

i 

£y& = Q. ft 
4e = 0.233 ft 

Dynamic pressure, air density and air density ratio are as follows: 

q= 0.5PV 2 

p = (o)(.0023769) slugs/ft 3 
a = Standard air density ratio at current altitude 
p 0 = (,629)(.0023769) slugs/ft 3 (i.e., h=15,000 ft) 

The gravity constant g is 32.174 ft/sec 2 - 
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2.2 Angle of Attack, Sideslip and Total Speed 

With respect to body axes, the angle of attack, a, the sideslips angle, P, and the 
total speed, V, are defined as 


a=mn '{^) 


P=Sin l \- 


V 2 = u 2 + v 2 + w 2 


- n<a<n 


-n<P<n 


(26) 

(27) 

(28) 


Our high alpha research vehicle (HARV) analytical model is derived from the wind- 
tunnel model described in 1 1 1. It is a full, nonlinear 6 DOF, rigid-body dynamic model 
whose aerodynamic forces and moments are calculated from a large wind-tunnel-derived 
data base using table look-ups with linear interpolation. The range of angle of attack, 
sideslip angle, mach number and altitude in that model are given in Table 2.2. 


Table 2.2. 

Range of State Variables in Aerodynamic Coefficients 

a 

-10° to 90« 

angle of attack (alpha) 

P 

-20° to 20° 

sideslip (beta) 

I M 

5 

0.2 to 2.0 

Mach number 

h 

0 to 60,000 ft 

altitude 

The aerodynamic coefficients are functions of the following control variables as 
well as angle of attack, sideslip, Mach number, altitude, roll, pitch and yaw rates: Aileron 


deflection (5a), Rudder deflection (5r) and Stabilator deflection (8h ). Throttle (5j) is an 
engine control variable. These control variables and their limits are given in Table 2.3. 
The effects of leading edge flap, trailing edge flap, speed brake, landing gear, etc are not 
considered. 


Table 2.3 Control Variables and Their Limits 

8a Aileron deflection (-25°> 25°) 

8r Rudder deflection (-30°, 30°) 

8h Stabilator deflection (-24°, 10.5°) 

8t Throttle 30° (idle) to 131° (full after burner) 
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2.3 Mathematical Structure of Aerodynamic Coefficients 

We consider the following state and control dependency structure for the 
coefficients of the aerodynamic model of the wind-tunnel based high angle of attack vehicle 
model. The effects of leading edge flap, trailing edge flap, speed brake, landing gear, etc 
are not considered. Only lift and pitching moment coefficients have an unsteady flow 
portion; they include the effect of the time rate-of-change of angle of attack. The other 
coefficients only have a steady flow part; they are explicit functions of airplane velocity 
states and control surface positions. 

Drag ('d 

C D =C D (a, M,h, 5h) (29) 

Lift Cl 


C L = C Lo (a,M,h,Sh) + ^[C Lf ( cc,M,h,)q + C Ld ( a,M,h)d\ 


(30) 


Pitching Moment C m 

C n = C mg (a, M,h, 8h ) + ^[c m? (a,M,% + C m ^ ( a,M,h)oc 
Side Force C y 

C =C (a,p,M,h,8a,8r) + C y (oc,M,h)P + 

y * n r 


(31) 


b_ 

2V 



(a,M,h)p + 


C y «x,M,h) 


r 


(32) 


Moment C / 


C^ = C^ (a,p,M,h,8a,Sr) + C^ (a,M,h)p 

° P 


+- 


2V 


(a,M,h)p + C^ (a,M,h) r 
. P r 


(33) 
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Yawing Moment C n 


C = C (a,p,M,h,8a,8r,8h) + C (a,M,h)f5 + 
nn o n p 


b_ 

2V 


C (a,M,h)p + C n (a,M,h) r 
L P r 


( 34 ) 


2.4 Explicit Equations of Motion 

Since the lift and pitching moment coefficients depend on a the right-hand sides of 
Eqs. (1), (3), (4)-(6) depend on the differentials u and vv through the relation 

uw - iiw 
a = — T 

u + w ( 35 ) 

which follows from the definition of a. 


Substitution of expressions for X, D, L, C D and C L into the u Eq.(l) yields 


F +C,sin(a)uw 

— h - J 


1 + C^sin(a)w 


in which 


F = rv-qw- gsin(d)- F D cos(a) + — 


m 


qS . 

+— sin(a) 
m 


C. +— C, q 

L 2V L 
o q 


(36) 


(37) 
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(38) 


D 


VSC D 

(i 

m 


h m(u 2 +w 2 ) 


Substitution of expressions for Z, D, L . Co, Cl and 


C, cos (a)w F 

F +— - 

w 1 + C^ sin(a)w 


B 

R 


in which 

T 

F = qu- pv + g cos( d) cos(O)- F sin( a) + — 

W L) m 

qS 

— - — cos(a) 
m 


C, +— C. q 
L 2VB 
° R J 


B 

q 


= 1 + C^cos(a)u- 


cos( a )wC^ sin( a)u 
1 + C^sin(a)w 


(39) 

i into the w Eq. (3) yields 

(40) 


(41) 

(42a) 
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or, equivalently. 


C, cos(a)u 

B =1+ — 

<7 1 + sin( a )w 


(42b) 


Substitution of w, Eq. (40), into u , Eq. (36), yields the right-hand-side of u 
withoui differentials. 

Since u and w have been defined so that their right-hand sides are free of 
differei tials, these expressions, Eqs. (36) and (40), can be substituted into Eq. (35) to give 
a free of differentials in its right-hand side. Consequently, expressions (35), (36) and 

(40) co istitute explicit expressions for a, u, and vv which are free of differentials in their 
right-h nd sides. 


3. Analytical Model for Drag Coefficient 


The analytical model of the drag coefficient Cq is taken from the wind-tunnel model 
at an altitude h= 15,000 feet and a Mach number M=0.6. The analytical model for Co 0 is 

constructed at stabilator deflections 5h = 10.5°, 0°, -5° and -24°. The analytical models 
are functions of a from 0° to 90°; they are defined in Tables 3.1 and 3.2. The analytical 
formulae are presented in Table 3.3. Comparisons of the analytical models with the 
corresponding wind-tunnel model data are shown in Figure 3.1. The analytical models 
are also given in the computer code listing contained in Appendix B. 

Table 3.1 Definitions of Cd Analytical Models 


CD 0 (a, M=0.6, 8h= 10.5°, h=15,000 ft ) 
Cd 0 («> M=0.6, 8h= 0°, h=15,000 ft ) 
CD 0 (a, M=0.6, 8h= -5°, h=15,000 ft ) 
CD 0 «*> M=0.6, 8h= -24°, h=15,000 ft ) 


= CD0X(a) 

= CD0Z(a) 

= CD0N5(a) 
= CD0N(a) 


Table 3.2 Drag Coefficient Analytical Models 


Sh 

Mach Number = 0.6 

Sh = 10.5 

CD0X(a) 

o 

ii 

£ 

CD0Z(a) 

8h= -5 

CD0N5(a) 

8h = -24 

CD0N(a) 


1 2 








Table 3.3 Formulas for CD Model 


CD0X(a°) = (.4/2.75)tan-i(((78/80)a°+7)l/30) 
+(.6/2.75)tair 1 (-((78/80)a°+2)l/8) 
+(-.3/2.75)tair , (-((78/80)a°+5)l/90) 
+(-.2/2.75)tan-'(-((78/80)a°-6)l/5) 
+(1.95/2.75)tan- 1 (((78/80)a°-28)l/15) 

+(2.2/2. 75)tair 1 (((78/80)a°-58)l/40) 

+( 1.4/2. 75)tair 1 (-((78/80)a°-73) 1/30) 
+(2.3/2.75)tan* , (-((78/80)a°- 138) 1/20)-. 147 

CDOZ(aO) = (2.17/2. 10)[(.6/2.75)tair 1 (((77/80)a°+6) 1/30) 
+(.6/2.75)tan* 1 (-((77/80)a°+ 1)1/8) 

+(-. 3/2. 75)tan* 1 (-((77/80)a°+4) 1/90) 
+(-.2/2.75)tan- , (-((77/80)a°-7)l/10) 
+(1.95/2.75)tan-i|((77/80)a°-29)l/15) 

+(2.2/2. 75)tair 1 (((77/80)a°-59)l/40) 

+(1. 55/2. 75)tair 1 (-((77/80)a°-74) 1/30) 

+(2.3/2.75)tan 1 (-((77/80)a°-139)l/20)-.2834] + .0199 

CD0N5(a°) = (.32/2.75)tan- 1 (((80/85)a°+8)l/30) 

+(.6/2.75)tan-»(-( (80/85)cc°+3.5)l/6.5) 
+(-.3/2.75)tan' , (-l(80/85)a°+7)l/90) 
+(-.2/2.75)tan* 1 (-«(80/85)a°-4)l/15) 
+(1.95/2.75)tair I (((80/85)a°-28)l/15) 

+(2.25/2. 75) tan* 1 * ((80/85)a°-68)l/40) 

+ (1.664/2.75)tair l (-((80/85)a°-90)l/30) 
+(2.35/2.75)tair 1 i-((80/85)a°-140) 1/20)-. 246 j 

CD0N(aO) = (.s/2.75)tan- , ((a°+5)l/30) 

+(.6/2.75)tair 1 (-(u°-0)l/6) 

+(-.25/2.75)tan- 1 (-(a°+3)l/90) 

+(-. 15/2. 75)tair 1 (-(a°-4) 1/40) 
+(1.85/2.75)tan* 1 ((a°-30)l/28) 

+(2. 3/2.75) tan’ 1 ((a o -60) 1/40) 
+(1.15/2.75)tan- 1 (-(a°-85)l/30) 

+(2.3/2.75)tair 1 (-(a°- 140) 1/20)-. 2425 
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CDO • C QQ 


Analytical Mcxlel of C D 

Wind-Tunnel Analytical 


CD OX (a) 


CDOZ(a) 










4. Analytical Model for Lift Coefficient 


The analytical model of the lift coefficient Cl is taken from the wind-tunnel model 
at an altitude h=15,000 feet and Mach numbers M=0.6 and 0.9. The analytical model for 
C Lo is constructed at stabilator deflections 8h = 10.5° and -24°. The analytical models 

are functions of a from 0° to 90°; they are defined in Tables 4.1 and 4.2. The analytical 
models of pitch and alpha dot rates are constructed at 0.6 Mach. The analytical formulae 
are presented in Table 4.3. Comparisons of the analytical models with the corresponding 
wind-tunnel model data are shown in Figures 4.1 and 4.2. The analytical models are also 
given in the computer code listing contained in Appendix B. 

Table 4.1 Definitions of Cl Analytical Models 


CLo(«> M=0.6, 8h= 10.5°, h=15,000 ft ) 
CL 0 («> M=0.9, 8h= 10.5°, h=15,000 ft ) 
CL 0 («> M=0.6, 8h= -24°, h=15,000 ft ) 
CL 0 (<x, M=0.9, 8h= -24°, h=15,000 ft ) 
CL q (a, M=0.6, h=15,000 ft ) 
CL adot (a, M=0 * 6 ’ ^=15,000 ft ) 


= CL0X6(a) 
= CL0X9(a) 
= CL0N6(a) 
= CL0N9(a) 
= CLQ(a) 

= CLAD(a) 


Table 4.2 Lift Coefficient Analytical Models 


Sh 

Mach Number 

0.6 

0.9 

8h= 10.5 

CL0X6(a) 

CL0X9(a) 

Sh = -24 

CL0N6(a) 

CL0N9(a) 


1 5 










Table 4.3 Formulas for Cl Model 

CL0X6(a°) = (.86/2.75)tan»(-(a°+5) 1/100) 
+(2.19/2.75)tair 1 ((a°-5)l/7) 
+(.9/2.75)tair 1 ((a°-24)l/17) 
+(1.71/2.75)tan' , (-(a°-53)2/25) 
+(.41/2.75)tair 1 (-(a°-70)2/7)-.95 

CL0X9(a°) = (.86/2.75) tair 1 (-(a°+5) 1/100) 
+(2.59/2.75)tan- 1 ((aO-3)l/7) 
+(1.6/2.75)tanr I ((a°-20)l/22) 
+(3.41/2.75)tair 1 (-(a°-57)l/30) 
+(.41/2.75)tan- 1 (-(a°-70)l/20)-.65 

CL0N6(a°) = (1.06/2.75)tan* 1 (-(a°+5)l/100) 

+( 1. 79/2. 75)tan- 1 ((a°-5) 1/7) 
+(2.5/2.75)tan' 1 ((a°-15)l/22) 

+(2.7 l/2.75)tair I (-(a°-59) 1/50) 
+(1.21/2.75)tan- , (-(a°-70)l/20)-.72 

CL0N9(a°) = (1.06/2.75)tan‘ 1 (-(a°+5)l/100) 
+(1.79/2.75)tan-i((a°-5)l/7) 
+(2.5/2.75)tan* 1 ((a°-13)l/22) 
+(2.71/2.75)tan' 1 (-(a°-59)l/50) 

+( 1. 21/2. 75)tair 1 (-(a°-70) 1/20)-. 80 

CLQ(a°) = (.26/2.75)tan* 1 (-(a o -5)l/10) 
+(-2.39/2. 75)tan- 1 ((aO-6)l/3) 
+(2.4/2.75)tan- 1 ((a°-15.5)l/3) 

+ (2.0/2.75)tan-»(-(a°-20)l/5) 
+(4.3/2.75)tair 1 ((a°-37)l/4.5) 

+(2. 2/2. 75)tair 1 (-(a°-45) 1/15) 
+(2.2/2.75)tan* 1 (-(a°-80)l/15) 
+(-.45/2.75)tan' , (-(a°-76)l/3.5)+4.2 

CLAD(a°) = (1.32/jc)tair , (-(a°-5)57c/18) 
+(-.75)tan-i((a°-45)l/2)+1.8 
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CLO 


Analytical Model of C L 


Wind-Tunnel 


Analytical 


CL0X6(a) CL0X9(a) 




8h= 10.5°, M=0.6 


8h= 10.5°,M=0.9 


CL0N6(a) 


CL0N9(a) 




5h= -240, m=0.6 


8h= -24°,M=0.9 


Figure 4.1: Comparison of Wind-Tunnel and Analytical Lift Coefficient Cj^ for 
h=15,000 feet: 8h=10.5°, -24° and M=0.6,0.9. 
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Analytical Model of C L 


Wind-Tunnel 


Analytical 


CLQ(a) 


CLAD(a) 



h= 15,000 feet, M=0.6 h= 15,000 feet, M=0.6 



Figure 4.2: Comparison of Wind-Tunnel and Analytical Lift Coefficient Derivatives C Lc 
and C Lad ^ fcr h=15, 000 feet and M=0.6. 
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5 . Analytical Model for Pitching Moment Coefficient 


The analytical model of the pitching moment coefficient C m is taken from the wind- 
tunnel model at an altitude h=15,000 feet. The analytical model for C mQ is constructed at 

stabilator deflections 8h = 10.5°, 5° , 2°, ()«, -5°, -12.5° and -24° for the Mach numbers 
M=0.3, 0.6, 0.8 and 0.9. The analytical models of pitch and alpha dot rates are 
constructed at 0.6 Mach. The analytical models are functions of a from 0° to 90°; they 
are defined in Tables 5.1a, b,c,d and 5.2. The analytical formulae are presented in Tables 
5.3 a,b,c,d. Comparisons of the analytical models with the corresponding wind-tunnel 
model data are shown in Figures 5.1 to 5.8. The pitch and alpha dot rates derivatives 
C mq arK * ^m a dot 3X0 given in Figure 5.8. The analytical models are also given in the 
computer code listing contained in Appendix B . 


Table 5.1a Definitions of C m Analytical Models at 0.6 Mach 


C mo («> M=0.6, 8h= 10.5°, h=15,000 ft ) 
Cm o (a,M=0.6, 5h= 5°, h=15,000 ft ) 
C mo (a. M=0.6, 5h= 2°, h=15,000 ft ) 
Cm 0 (a, M=0.6, 5h= 0°, h=15,000 ft ) 
Cm 0 (a> M=0.6, 8h= -5°, h=15,000 ft ) 
Cmo(o, M=0.6, 8h= -12.5°, h=15,000 ft ) 
C mo (a, M=0.6, 8h= -24°, h=15,000 ft ) 
C mq (a, M=0.6, h=15,000 ft ) 
c m a dot (M=0.6, h=15,000 ft ) 


= CmOX6(a) 

= CmOX56(a) 
= CmOX26(a) 
= Cm006(a) 

= CmON56(a) 
= CmOZ6(a) 
= CmON6(a) 
= CMQ(a) 

= CMAD(a) 


Table 5.1b Definitions of C m Analytical Models at 0.8 Mach 


C mo (a» M=0.8, 8h= 10.5°, h=15,000 ft ) 
C mo (a,M=0.8, 8h= 5°, h=15,000 ft ) 
C mo (a, M=0.8, 8h= 2°, h=15,000 ft ) 
C mo (a, M=0.8, 8h= 0°, h=15,000 ft ) 
C mo (a, M=0.8, 8h= -5°, h=15,000 ft ) 
C mo (a, M=0.8, 8h= -12.5°, h=15,000 ft ) 
C mo (a, M=0.8, 8h= -24°, h=15,000 ft ) 


= CmOX8(a) 

= CmOX58(a) 
= CmOX28(a) 
= Cm0X08(a) 
= CmON58(a) 
= CmONZ8(a) 
= CmON8(a) 


Table 5.1c Definitions of C m Analytical Models at 0.9 Mach 


C mo (a, M=0.9, 8h= 10.5°, h=15,000 ft ) 
Cm 0 (a,M=0.9, Sh= 5°, h=15,000 ft ) 
C mo (a, M=0.9, 8h= 2°, h=15,000 ft ) 
C mo (a, M=0.9, 8h= 0°, h=15,000 ft ) 
C mo (a, M=0.9, 8h= -5°, h=15,000 ft ) 
C mo (a, M=0.9, 8h= -12.5°, h=15,000 ft ) 
C mo (a, M=0.9, 8h= -24°, h=15,000 ft ) 


= CmOX9(a) 

= CmOX59(a) 
= CmOX29(a) 
= Cm0X09(a) 
= CmON59(a) 
= CmONZ9(a) 
= CmON9(a) 
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Table 5.1d Definitions of Cm Analytical Models at 0.3 Mach 


C mo (a, M=0.3, 8h= 10.5°, h=15,000 ft ) = CmOX3(a) 

C mo (a, M=0.3, 8h= -12.5°, h=15,000 ft ) = CtnONZ3(a) 
C ra# (a, M=0.3, 8h= -24°, h=15,000 ft ) = CmON3(a) 

Table 5.2 Pitching Moment Coefficient Analytical Models 


Mach Number 


Sh 



0.3 

0.6 

0.8 

0.9 

10.5 

CmOX3(a) 

CmOX6(a) 

CmOX8(a) 

CmOX9(a) 

5 

CmOX56(a) 

CmOX56(a) 

CmOX58(a) 

CmOX59(a) 

2 

CmOX26(a) 

CmOX26(a) 

CmOX28(a) 

CmOX29(a) ) 

0 

Cm006( a) 

Cm006( a) 

Cm0X08(a) 

Cm0X09(a) 

-5 

CmON56(a) 

CmON56(a) 

CmON58(a) 

CmON59(a) 

1 

CmONZ3(a) 

CmOZ6(a) 

CmONZ8(a) 

CmONZ9(a) 

-24 

CmON3(a) 

CmON6(a) 

CmON8(a) 

CmON9(a) 
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Table 5.3a Formulas for Cm Model at 0.6 Mach 


CmOX6(a°) = (.26/2.75)tan- 1 (-(ao.5)l/10)+(-J9/2.75)tair 1 ((a°-l)l/8) 
+(.8/2.75)tan- 1 ((a°-5)l/13)+(.70/2.75)tan- 1 (-(a°-10)l/65) 
+(1.2/2.75)tan- , ((a 0 -49)l/15)+(2.1/2.75)tan' 1 (-(a°-69)l/15) 
+(-.45/2.75)tair , (-(a°-77)l/2)-.398 

CmOX56(a°) = (.26/2.75)tan- 1 (-(a o -5)l/10)+(-.39/2.75)tan 1 ((a°-l)l/12) 
+(.9/2. 75)tan- 1 ((a°-5) l/15)+(. 85/2. 75)tan- 1 (-(a°- 10) 1/30) 
+(1.35/2.75)tan' 1 ((a°-49)l/19)+(2.2/2.75)tan* 1 (-(a°-69)l/15) 
+(-. 45/2. 75)tan' 1 (-(a°-77) 1/2)-. 368 

CmOX26(a°) = (.26/2.75)tan- 1 (-(a o -2)l/10)+(-.39/2.75)tan* 1 ((a°-l)l/10) 
+(1/2. 75)tair 1 ((a°-3)l/ll)+(. 85/2. 75)tair 1 (-(a°-7) 1/20) 
+(1.35/2. 75)tair 1 ((a°-51) l/19)+(2. 2/2. 75)tair 1 (-(a°-69) 1/1 5) 
+(-.45/2.75)tan- 1 (-(a°-77)l/2)-.31 

Cm006(a°) = (.36/2.75)tan' 1 (-(a°-5)l/30)+(-.29/2.75)tan' 1 ((a°-l)l/15) 
+(.9/2. 75)tan* 1 ((a°-5) l/35)+(. 80/2.75) tan‘ 1 (-(a°-48) 1/75) 
+(.9/2.75)tair , ((a°-52)l/10)+(2.1/2.75)tan- 1 (-(a°-69)l/15) 
+(-.45/2.75)tan->(-(a°-77)l/2)-.457 

CmON56(a°) = (.26/2.75)tan- 1 (-(a°-5)l/30)+(-.39/2.75)tan- 1 ((a°-l)l/30) 

+( 1.2/2. 75)tan- , ((a°-5) l/40)+(. 60/2. 75)tair , (-(oc°-8) 1/23) 

+( 1.3/2. 75)tan‘ 1 (-(ct°-60)l/65)+(2. 8/2. 75)tan' 1 ((a°-72) 1/55) 
+(2.3/2. 75)tair 1 (-(a°-73)l/19)+(-. 45/2. 75)tair , (-(a°-77) 1/2) 
-.188 

CmOZ6(a°) = (.26/2.75)tan-H-(a°-5)l/60)+(-.39/2.75)tan- 1 ((a o -l)l/14) 

+(. 8/2. 75)tan -1 ((a 0 -5)l/42)+(. 80/2. 75)tan- 1 (-(a°-20)l/55) 

+( 1.8/2. 75)tan* 1 ((a°-65)l/60)+(2. 4/2. 75)tair 1 (-(a°-69) 1/20) 
+(-.45/2.75)tan* 1 (-(ot°-79)l/2)-.158 

CmON6(a°) = (.26/2.75)tan-i(-(a°-5)l/60)+(-.39/2.75)tan- 1 ((a°-l)l/30) 
+(.8/2.75)tan- , ((a°-5)l/45)+(.80/2.75)tan* 1 (-(a°-l 0)1/65) 
+(1.8/2.75)tan- , ((a°-51)l/45)+(2.8/2.75)tan- , (-(a°-69)l/23) 
+(-.45/2.75)tair , (-(a°-79)l/2)-.148 

CMQ(a°) = (-.82/7i)tan- , (-(a°-5)2jc/18)+(2)tan- 1 (-(a°-32)l/6) 

+(4.55)tan-i((a°-43)3.5/l)+(-3.5)tan- , ((a 0 -57)l/5)-5.8 

CM AD(a°) = (-. 02 /jc)tan‘ , (-(a°-l) 57 c/ 18 )+(. 5 )tan- , ((a°- 6 ) 5 / 1 ) 
+(-.8)tan-i((a°-18)l/2)+(0.9)tan- 1 ((a°-45)l/2)-.9 
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Table 5.3b Formulas for Cm Model at 0.8 Mach 

CmOX8(a°) = (.26/2.75)tan- 1 ("(aO-5)l/15)+(-.33/2.75)tan- 1 ((a°-l)l/7) 

+(. 72/2. 75)tan- 1 ((a°-5)l/15)+(. 70/2. 75)tair 1 (-(a°-35) 1/75) 

+(1. 13/2. 75)tair 1 ((a°-51) 1/1 l)+(2.08/2.75)tair 1 (-(a°-67)l/17) 
+(-.45/2.75)tan-i(-(a°-78)l/4)-.440 

CmOX58(a°) = (.26/2.75)tan-i(-(a°-5)l/10)+(-.39/2.75)tan*i((a°-l)l/12) 

+(. 7/2. 75)tair 1 ((a°-5)l/15)+(. 85/2. 75)tan‘ 1 (-(a°-20) 1/40) 
+(1.45/2. 75)tan- 1 ((a°-52)l/15)+(2. 2/2. 75)tan* 1 (“(a°-69) 1/15) 
+(-.45/2.75)tan- , (-(a°-77)l/3.4)-.338 

CmOX28(a°) = (.36/2.75)tan-i(-(a°-5)l/35)+(-.29/2.75)tan-i((a°-l)l/30) 
+(l/2.75)tan- 1 ((a°-15)l/90)+(.75/2.75)tan' 1 (-(«°-48)l/110) 
+(.9/2.75)tan-M(a°-52)l/9)+(2.1/2.75)tan-i(-(a°-69)l/17) 
+(-.45/2.75)tan* 1 (-(a°-77)l/3)-.387 

Cm0X08(a°) = (.36/2.75)tan-i(-(a°-5)l/35)+(-.29/2.75)taiH((a°-l)l/30) 
+(l/2.75)tair 1 ((a°-15)l/90)+(.80/2.75)tan-i(-(a°-48)l/90) 

+ (.9/2.75)tan’ 1 ((tx°-52)l/9)+(2.1/2.75)tan- 1 (-(a°-69)l/17) 

+(-. 45/2. 75)tair 1 (-(a°-77) 1/3)-. 387 

CmON58(a°) = (.36/2.75)tan-i(-(a°-5)l/35)+(-.29/2.75)tan-i((a<>-l)l/30) 
+(1.3/2.75)tan- 1 ((a°-15)l/95)+(.80/2.75)tan- 1 (-(a°-47)l/35) 
+(l/2.75)tan- 1 ((a°-53)l/10)+(2.15/2.75)tan-»(-(a°-69)l/18) 
+(-.45/2.75)tair 1 (-(a°-78)l/2)-.367 

CmONZ8(a°) = (.36/2.75)tairi(-(aO-5)l/35)+(-.29/2.75)tan»((ao.l)l/30) 

+(1.25/2.75)tan-i((aO-15)l/95)+(.80/2.75)tan-i(-(a°-47)l/28) 
+ (l/2.75)tan-i((a°-53)l/10)+(2.25/2.75)tan-i(-(a°- 69)1/18) 
+(-.45/2.75)tan-i(-(a°-78)l/2)-.317 

CmON8(a°) = (.26/2.75)tanU-(a o -5)l/40)+(-.45/2.75)tan , ((a o -4)l/30) 
+(.7/2.75)tan , ((oc°-2)l/40)+(.80/2.75)tan- l (-(a°-37)l/25) 

+( 1.9/2.75)tarr 1 ((a°-52) l/25)+(2. 7/2. 75)tan* 1 (-(a 0 -69) 1 / 20 ) 

+ (-.45/2.75)tan-i(-(aO-79)l/2)-.198 
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Table 5.3c Formulas for C m Model at 0.9 Mach 

CmOX9(«°) = (.26/2.75)tan- 1 (-(a°-5)l/10)+(-.39/2.75)tan- 1 ((a°-0)l/10) 
+(l/2.75)tan-i((a°-3)l/25)+(.70/2.75)tan-i(-(a°-7)l/25) 

+(1. 3/2. 75)tair 1 ((a°-50) 1/1 6)+(2. 1/2. 75)tair , (-(a°-69) 1/15) 
+(-.45/2. 75)tan->(-(a°- 76)1 /3.5)-.433 

CmOX59(a°) = (.26/2.75)tan-i(-(a°-5)l/4)+(-.39/2.75)tair 1 ((ao+2)l/5) 
+(l/2.75)tair 1 ((a°-l)l/15)+(.60/2.75)tan' 1 (-(oc°-18)l/13) 
+(1.3/2.75)tan-i((aO-51)l/14)+(2.15/2.75)tan-i(-(a°-69)l/15) 
+(-.45/2.75)tan-i(-(a°-76)l/3.5)-.490 

CmOX29(a°) = (.26/2.75)tan* 1 (-(a°-5)l/10)+(-.39/2.75)tair 1 ((a°+4)l/7) 
+(1.33/2. 75)tan* 1 ((a°-l)l/40)+(. 60/2. 75)tair 1 (-(a°- 19) 1/17) 
+(1.9/2.75)tair 1 ((a°-51)l/9)+(2.05/2.75)tair 1 (-(a°-69)l/15) 
+(-.45/2.75)tair 1 (-(a°-76)l/3.5)-.450 

Cm0X09(a°) = (.26/2.75)tair 1 (-(a°-5)l/10)+(-.39/2.75)tan-i((a°+2)l/10) 
+(1.3/2.75)tan- , ((a°-l)l/40)+(.60/2.75)tan- 1 (-(a°-19)l/15) 
+(.9/2.75)tan- I ((a°-51)l/9)+(2.11/2.75)tan- , (-(a°-69)l/15) 
+(-.45/2.75)tan- 1 (-(a°-76)l/3.5)-.490 

CmON59(a°) = (.26/2.75)tan-»(-(a<>-5)l/ll)+(-.39/2.75)tan-i((ao + 2)l/10) 
+(1.5/2.75)tair 1 ((a°-l)l.l/48)+(.60/2.75)tair 1 (-(a°-19)l/15) 
+(. 8/2. 75)tair 1 ((a°-51)l/10)+(2. 32/2. 75)tair 1 (-(a°-70) 1 / 20 ) 
+(-.45/2.75)tan" 1 (-(<x°-78)l/3.5)-.490 

CmONZ9(a°) = (.26/2.75)tan- 1 (-(a°-5)l/7)+(-.39/2.75)tan- 1 ((a°-0)l/10) 

+(1. 6/2. 75)tan* 1 ((a°-5)l/50)+(. 60/2. 75) tair 1 (-(a°-32) 1/11) 

+(.8/2.75)tan-i((aO-51)l/ll)+(2.23/2.75)tan-i(-(aO-69)l/19) 

+(-.45/2.75)tair 1 (-(a°-78)l/3.5)-.410 

CmON9(a°) = (.16/2.75)tan- 1 (-(a°-5)l/40)+(-.39/2.75)tan I ((a°-3)l/8) 

+(1.2/2.75)tan* , ((a°-15)l/120)+(.70/2.75)tair 1 (-(a°-25)l/50) 

+(2/2.75)tan' 1 ((a°-52)l/75)+(2/2.75)tair , (-(a°-69)l/20) 

+(-. 42/2. 75)tan _1 (-(a°-79) 1/4)-. 068 
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Table 5.3d Formulas for C m Model at 0.3 Mach 


CmOX3(a°) = (.26/2.75)tan* 1 (-(a o -5)l/10)+(-.39/2.75)tair 1 ((a o - 1)1/8) 

+(. 75/2. 75)tan- 1 ((a°-5)l/13)+(. 70/2. 75)tan- 1 (-(a° -10) 1/65) 

+(1.2/2.75)tan* 1 ((a°-49)l/15)+(2.1/2.75)tan‘ 1 (-(a°-69)l/15) 

+(-.45/2.75)tan* 1 (-(a°-77)l/2)-.398 

CmONZ3(a°) = (.26/2.75)tan-i(-(aO-5)l/60)+(-.39/2.75)tan-i((a°-l)l/14) 
+(. 85/2. 75)tan' 1 ((a 0 -5)l/42)+(. 80/2. 75)tan* 1 (-(a°-50) 1 / 60 ) 
+(1.8/2.75)tair I ((a o -70)l/54)+(2.4/2.75)tan* 1 (-(a°-69)l/25) 
+(-. 45/2. 75)tan- , (-(a°-79) 1/2)-. 158 

CmON3(a°) = (.26/2.75)tarr 1 (-(a°-5) l/60)+(-. 39/2. 75)tan' 1 ((a°- 1)1/30) 
+(. 8/2. 75)tair 1 ((a°-5)l/45)+(. 80/2. 75)tair 1 (-(cc°- 10) 1/65) 

+ ( 1.8/2. 75)tan’ 1 ((a°-49)l/40)+ (2. 8/2. 75)tan- 1 (-(a°-69) 1/23) 
+(-.45/2. 75)tan- 1 (-(a°-79)l/2)-. 138 

CmON53(a°) = CmON56(a°) 

Cm0X03(a°) = Cm006(a°) 

Cm0X23(a°) = CmOX26(a°) 

CmOX53(a°) = CmOX56(a°) 
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Figure 5.1: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C 
for h=15,000 feet and 5h=10.5°. 
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Analytical Model of C m 


Wind-Tunnel 


Analytical 


CM0X56(a) 


CM0X56(a) 




M=0.3 


M=0.6 


CM0X58(a) CM0X59(a) 




M=0.8 


M=0.9 


Figure 5.2: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C m ^ 
for h= 15,000 feet and 5h=5°. 
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Analytical Model of C m 

Wind-Tunnel Analytical 


CM0X26(a) CM0X26(a) 




M=0.3 


M=0.6 


CM0X28(a) 


CM0X29(a) 




M=0.8 M=0.9 


Figure 5.3: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C m ^ 
for h= 15,000 feet and 5h=2°. 
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Analytical Model of C m 
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Analytical 


CM006(a) 


CM006(a) 




M=0.3 


M=0.6 


CM0X08(a) 


CM0X09(a) 




M=0.8 


M=0.9 


Figure 5.4: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C m ^ 
for h=15,000 feet and 8h=0°. 
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Wind-Tunnel Analytical 
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M=0.3 M=0.6 


CM0N58(a) (j:M0N59(a) 



M=0.8 M=0.9 


Figure 5.5: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C mQ 
for h= 15,000 feet and 8h=-5°. 
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Analytical Model of C m 
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M=0.3 


M=0.6 


CM0NZ8(a) 


CM0NZ9(a) 




M=0.8 


M=0.9 


Figure 5.6: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C mf 
for h=15,000 feet and 8h=-12.5°. 
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CM0N8(a) 


CM0N9(a) 




M=0.8 
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Figure 5.7: Comparison of Wind-Tunnel and Analytical Pitching Moment Coefficient C mQ 
for h=15,000 feet and 8h=-24°. 
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Analytical Model of C m 
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Figure 5.8: Comparison of Wind-Tunnel and Analytical Lift Coefficient Derivatives C n 

and C m , for h= 15,000 feet and M=0.6. 
m adot 
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6. Analytical Model for Side Force Coefficient 


The analytical model of the side force coefficient Cy is taken from the wind-tunnel 
model at an altitude h= 15,000 feet and a Mach number M=0.6. The analytical model for 
C yQ is constructed at P = 0°, 20°; 5a=+ 25°; and 8r = + 30°. The analytical models are 

functions of a from 0° to 90°; they are defined in Tables 6.1 and 6.2. The analytical 
formulae are presented in Table 6.3. Comparisons of the analytical models with the 
corresponding wind-tunnel model data are shown in Figures 6.1 to 6.3. The sideslip 
derivative C y p is taken as the constant .000206. The roll and yaw rate derivatives C yp 

and C yr are given in Figure 6.3. The analytical models are also given in the computer 
code listing contained in Appendix C. 


Table 6.1 Definitions of Cy Analytical Models 


C yo (a, p=20°, M=0.6, 8a= 25°, 8r=-30°, h=15,000 ft ) 
C yo (a, P=20°, M=0.6, 8a= 25°, 8r= 30°, h=15,000 ft ) 
C yo «x, (3=20°, M=0.6, 8a=-25°, 8r=-30°, h=15,000 ft ) 
C yo (a, 0=20°, M=0.6, 8a=-25°, 8r= 30°, h=15,000 ft ) 
C yo (a, p=0° , M=0.6, 8a= 25°, 8r=-30°, h=15,000 ft ) 
C yo (a, P=0° , M=0.6, Sa= 25°, 8r= 30°, h=15,000 ft ) 
Cy p (a, M=0.6, h=15,000 ft ) 

C yr (a, M=0.6, h=15,000 ft ) 

Cy p (M=0.6, h=15,000 ft ) 


= CyOXNB2(a) 
= CyOXXB2(a) 
= CyONNB2(a) 
= CyONXB2(a) 
= CyOXNBO(a) 
= CyOXXBO(a) 
= CYP(a) 

= CYR(a) 

= CYB(a) 


Table 6.2 Side Force Coefficient Analytical Models: M = 0.6 


Sideslip 

N. Rudder 

Aileron 

Sr = -30 

Sr = 30 

P = 20 

6u= 25 

CyOXNH2( a) 

CyOXXB2( a) 

& = -25 

CyONNB2( a) 

CyONXB2( a) 

p = o 

Su= 25 

CyOXNBOf a) 

Cy0XXB0( a) 

Sa = -25 

-Cy0XXB0( a) 

-CyOXNBO(o) 
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Table 6.3 Formulas for Cy Model 

CyOXNB2(a°) = (.H^jcHarr’ttaO-MHlx/lS) - .394 

CyOXXB2(a°) = (.06/n)tan-i(ao/3) 

+(.09/7C)tan- 1 ((a°-31)5/8) 

+ (.06/x)tan* 1 (-(a°-46)3/10) 
+(.03/ic)tan- 1 ((a°-63)3/7) 
+(.09/TC)tan- 1 (-(a°-75)3/5) 
+(.04/jc)tan- 1 ((a o -85)4/5) - .285 

CyONNB2(a°) = (.219544/7i)tan- 1 ((oc°-21)5/56) 
+(.0636375/jc)tan- , (-(a°-74)l/4) 

+ (.0709125/jc)tan-i((a°-85)3/10) - .36444 

CyONXB2(cc°) = (.047/x)tan-i(a°/2) 

+(.021/K)tan-i(-(a°-17)5/4) 

+ (.037/ji)tan-i((a°-32)5/3) 
+(.06/7t)tan-i(-(a°-76)5/4) 
+(.043/jc)tan- 1 (a°-85) - .248 

CyOXNBO(aO) = (.029624)tan-i((a°-20)4/25) 
+(-.0020868)tan-i(-(a°-12)5) 
+(6.99)exp(-(a°+17.6)) - .075535 

CyOXXBO(aO) = (.01216/2.75)tan-i(a°3/4) 

+(. 03247/2. 75)tair 1 (-(a°- 13) 1/4) 
+(.00891/2.75)tan' 1 ((a°-29. 5)2/3) 

+(. 03058/2. 75)tan' 1 (-(a°-46)2/5) 

+ (.02759/2.75) tan' 1 (-(a °-75)3/40) + .03477 

CyONXBO(a°) = - CyOXNBO(a°) 

CyONNBO(a°) = - CyOXXBO(a°) 

CYP(a°) = (.086/7i)tan' 1 (a°107c/18) 

+(.096/7C)tan- 1 (-(a°-23)107t/18) 

+ (.22/7t)tan* 1 (-(a°-45)107t/18) 
+(.256/7C)tan' 1 ((a°-54)107c/18) - .047 

CYR(a°) = (.17/7t)tan' 1 ((a°-4)107t/18) 
+(.55/7t)tan' 1 (-(a°-20)107c/18) 

+ (.54/7t)tan- 1 ((a°-45)107c/18) 
+(.26/7t)tan-t(-(a°-61)107c/18) 

+ .07 

CYB(a°) =.000206 
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CTO 


Analytical Model of Cy 


Wind-Tunnel 


Analytical 


CY0XNB2(a) CY0XXB2(a) 




8a= 25°, 8r= -30° 


8a= 25°, 8r= 30° 


CY0NNB2(a) CY0NXB2(a) 




8a= -25°, 8r= -30° 


8a= -25°, Sr= 30° 


Figure 6.1: Comparison of Wind-Tunnel and Analytical Side Force Coefficient C yQ for 
h= 15,000 feet: p=20° and M=0.6. 
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CYOXNBO(a) CYOXXBO(a) 




5a= 25°, 5r= -30° 5a= 25°, 6r= 30° 

Figure 6.2: Comparison of Wind-Tunnel and Analytical Side Force Coefficient C y ^ for 
h=15, 000 feet: p=0° and M=0.6. 

c y p (a) C yr (a) 




h= 15,000 feet, M=0.6 h= 15,000 feet, M=0.6 


Figure 6.3: Comparison of Wind-Tunnel and Analytical Side Force Coefficient Derivatives 
G. and C v for h=15, 000 feet and M=0.6. 

y p 
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7. Analytical Model for Rolling Moment Coefficient 


The analytical model of the rolling moment coefficient Q is taken from the wind- 
tunnel model at an altitude h=15,000 feet and Mach numbers M=0.6 and 0.9. The 
analytical model for Q 0 is constructed at P = 0°, 20°; 8a=+ 25°; and 8r = + 30°. The 

analytical models are functions of a from 0° to 90°; they are defined in Tables 7.1 a,b 
and 7.2 a,b. The analytical formulae are presented in Tables 7.3 a,b. Comparisons of 
the analytical models with the corresponding wind-tunnel model data are shown in 
Figures 7.1 to 7.4. The roll and yaw rate derivatives C/ p and C/ r and the sideslip 

derivative C ^ are given in Figure 7.4. The analytical models are also given in the 
computer code listing contained in Appendix C. 


Table 7.1a Definitions of C/ Analytical Models at 0.6 Mach 


C/ 0 (a, P = 20°, M=0.6, 8a= 25°, Sr=-30°, h=15,000 ft ) 
C/ 0 (a, P = 20°, M=0.6, 8a= 25°, 8r= 30°, h=15,000 ft ) 
C/ 0 (a, P = 20°, M=0.6, 8a=-25°, 8r=-30°, h=15,000 ft ) 
C/ 0 (a, P=20°, M=0.6, 8a=-25°, 8r= 30°, h=15,000 ft ) 
C/ 0 (a, P=0° , M=0.6, 8a= 25°, 8r=-30°, h=15,000 ft ) 
C/ 0 (a, P=0° , M=0.6, Sa= 25°, 8r= 30®, h=15,000 ft ) 
C/p (a, M=0.6, h=15,000 ft ) 

C/ r (a, M=0.6, h=15,000 ft ) 

C/p(M=0.6, h= 15,000 ft ) 


= C/0XNB2(a) 
= C/0XXB2(a) 
= C/0NNB2(a) 
= C/0NXB2(a) 
= C/0XNB0(a) 
= C/0XXB0(a) 
= LC/P(a) 

= C/R(a) 

= C/B(a) 


Table 7.1b Definitions of C/ Analytical Models at 0.9 Mach 


C/ 0 (a, P=20°, M=0.9, 8a= 25°, 8r=-30°, h=15,000 ft ) 

C/ 0 (a, P=20°, M=0.9, 8a= 25°, 8r= 30°, h=15,000 ft ) 

C/ 0 (a, P=20°, M=0.9, 8a=-25°, 8r=-30°, h=15,000 ft ) 
C/ 0 (a, P=20°, M=0.9, 8a=-25®, 8r= 30°, h=15,000 ft ) 
C/ 0 (a, P=0° , M=0.9, 8a= 25°, 8r=-30°, h=15,000 ft ) 

C/ 0 (a, p=0° , M=0.9, 8a= 25°, 8r= 30°, h=15,000 ft ) 

C/ p (a, M=0.9, h=15,000 ft ) 


= C/0XN2(a) 
= C/0XX2(a) 
= C/0NN2(a) 
= C/0NX2(a) 
= C/0XN0(a) 
= C/0XX0(a) 
= HQP(a) 



Tabic 7.2a Rolling Moment Coefficient Analytical Models: M = 0.6 


Sideslip 

N. Rudder 

Aileron 

8r = -30 

6r = 30 

P s 20 

&= 25 

CfDXNB2( a) 

C/0XXB2( a) 

6a = -25 

Cfl)NNB2( a) 

C«)NXB2( a) 

o 

ii 

ea 

6a= 25 

Ctf)XNB0( a) 

CA)XXB0( a) 

8a = -25 

-Cf0XXB0( a) 

-Cf0XNB0( a) 


Table 7.2b Rolling Moment Coefficient Analytical Models: M = 0.9 


Sideslip 

Rudder 

Aileron 

8r = -30 

6r = 30 

o 

CNJ 

II 

CO. 

&= 25 

C«XN2( a) 

C«XX2( a) 

6a = -25 

C20NN2( a) 

C/ONX2( a) 

p = o 

6u= 25 

C*)XN0( a) 

C«XX0( a) 

6a = -25 

-a0XX0( a) 

-Cf0XN0( a) 
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Table 7.3a Formulas for C / Model at 0.6 Mach 

C/0XNB2(a°) = (.OSS/Z^tan^-faQ-lS)^) - .0065 

C*0XXB2(a°) = (.02771/2.75)tan-i(-(a°-2)l/5) 

+(.06763/2. 75)tan' 1 (-(a°-19)l/4) 

+(. 006/2. 75)tan* 1 ((a°-22)2/l) + .00081 

C/0NNB2(a°) = (.0344/2.75)tan*»(-(a°-5)3/13) 

+(.037/2. 75) tan* , ((a°-24)2/5) 
+(.011/2.75)tan- 1 (-(a°-38)2/l) 

+(. 012/2. 75)tair 1 ((a°-42)7/4) 
+(.011/2.75)tair 1 (-(a°-52)3/8) 
+(.0176/2.75)tan-i((a°-73)3/13) - .0553 

C*0NXB2(a°) = (.0395/2.75)tan-i(-(a°-5)3/13) 

+(. 0295/2. 75) tan' 1 ((cx°-25)3/7) 
+(.0126/2.75)tan-!(-(a°-38)4/3) 
+(.0114/2.75)tair 1 ((a°-42)2/l) 
+(.0082/2.75)tan- 1 (-(a°-5 1)1/2) 

+(. 0132/2. 75)tan‘ 1 ((a°-70)2/5) - .0479 

C/0XNB0(a°) = (.0041/2.75)tan* 1 (-(a°-4)2/3) 
+(.003/2.75)tair 1 (-(a°-20)l/3) 
+(.011/2.75)tan- 1 (-(a°-59)2/5) 
+(-.00144/2.75)tan- 1 (-(a°+l)8/l) + .01793 

C*0XXB0(a°) = (.04226/2.75)tair 1 (-(a°-20)2/7) 
+(.00831/2. 75)tan' 1 (-(a°-53)4/7) 

+(. 00997/2. 75) tan' 1 ((oc°-65)4/5) 

+(.0101/2.75)tair 1 (-(a°-77.5)8/15) 

+(-.002/7C)tan 1 (-(a°-8)10)+.0286 


LCiP(a°) = (.15/7t)tan- , ((a°-12)107t/18) 

+(.25/7E)tan-i(-(a°-28)1007c/18) 

+(.55/7c)tan' 1 ((a°-41)100jc/18) 

+(.33/7i)tan' 1 (-(a°-50)107c/18)-.341 

C/R(a°) = (.304/ji)tan-i((a o -3)10jc/18) 
+(.22/7c)tair 1 (-(a°-50)2jc/18) 
+(-.026)exp((a°-85) 1/100)+. 018 

C/B(a°) = (.0001)[(6.32/7c)tan- , (-(aO-13)100ji/18) + 3.26] 
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Table 7.3b Formulas for Model at 0.9 Mach 

CiOXN2(a°) = (.048/2.75)tan-i(-(a°-18)8/92) 
+(.030/2.75)tair 1 (-(a°-23)8/92) 
+(.045/2.75)tan- 1 ((a°-80)8/92) - .0105 

C/0XX2(a°) = (.085/2.75)tan- 1 (-(a°-15)8/92) - .0198 

C/0NN2(a°) = -[(.0544/2.75)tan* 1 (-(a°-7)3/13) 

+(. 087/2. 75)tan‘ 1 ((a°- 17)2/7) 

+(. 008/2. 75)tair 1 (-(a°-25)2/l) 

+(.019/2.75)tan- 1 ((a°-28)7/4) 

+(.051/2.75)tan- 1 (-(a°-42)3/8) 

+(. 0096/2. 75)tan- 1 ((a°-55)3/13) - .05531(1/2.6) 

- .074 

C/0NX2(a°) = (.0295/2.75)tan- 1 (-(a°-25)3/13) 

+(.0295/2. 75)tan- 1 ((a°-42.5)3/7) 

+ (.0086/2. 75)tan- 1 (-(a°-15)4/3) 
+(.0014/2.75)tan' 1 ((ot°-42)2/l) 
+(.0162/2.75)tair 1 (-(a°-5 1)1/2) 

+ (.0012/2. 75)tan- 1 ((a°-70)2/5) - .0466 

C*0XN0(a°) = [(.0041/2.75)tair 1 ((ot°-8)2/3) 
+(.012/2.75)tan-i(-(a°-ll)l/3) 
+(.010/2.75)tair , (-(a°-15)l/3) 
+(.012/2.75)tair 1 (-(a°-35)l/3) 
+(.005/2.75)tan' 1 (-(a°- 100)2/5) 

+ (-.00144/2. 75)tarr 1 (-(a°-2)8) + .017931(1/2) 

- .0032 

C/0XX0(a°) = [(.04226/2.75)tair 1 (-(a°-8. 5)2/7) 

+(. 0103 1/2. 75)tarr 1 (-(a°-35)4/7) 

+(. 00997/2. 75)tair 1 ((a°-65)4/5) 

+ (.0131/2.75)tan-!(-(a o -77. 5)8/15) 
+(-.002/TC)tan- 1 (-(a°-8)10) + .0286](1/1.8) 


HC/P(a°) = (.28/7C)tan-i((a°-10)10jc/18) 

+(.25/ji)tair 1 (-(oc°-41)1007c/18) 
+(.55/7C)tan- 1 ((a°-41)1007t/18) 
+(.33/7c)tan- 1 (-(a°-50)107t/18) - .471 
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Analytical Model of C t 


Wind-Tunnel 
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C*0XNB2(a) CiOXXB2(oc) 




8a= 25°, 8r= -30° 


8a= 25°, 8r= 30° 


C£ONNB2(oc) 


CiONXB2(a) 




8a= -25°, 5r= -30° 


8a= -25°, 8r= 30° 


Figure 7.1: Comparison of Wind-Tunnel and Analytical Rolling Moment Coefficient C^ 
for h=15,000 feet , M=0.6 and p=20°. 
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Figure 7.2: Comparison of Wind-Tunnel and Analytical Rolling Moment Coefficient C 


e 0 


for h= 15,000 feet , M=0.9 and p=20°. 
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010 Oio 


Analytical Model of 


Wind-Tunnel 


Analytical 


CfOXNBO(a) 


CfOXXBO(a) 




8r= -30°, M=0.6 8r= 30°, M=0.6 


C^OXNO(a) 


aOXXO(a) 




8r= -30°, M=0.9 


8r= 30°, M=0.9 


Figure 7.3: Comparison of Wind-Tunnel and Analytical Rolling Moment Coefficient 
for h= 15,000 feet , Sa= 25° and p=0°. 



43 






HC1P 
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LCiP(a) [ C /p ] C^R(a) [ C t ] 
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M=0.6 
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Figure 7.4: Comparison of Wind-Tunnel and Analytical Rolling Moment Coefficient 
Derivatives C/ , C, andC/ for h= 15,000 feet andM=0.6, 0.9. 

<p < r tp 
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8. Analytical Model for Yawing Moment Coefficient 


The analytical model of the yawing moment coefficient C n is taken from the wind- 
tunnel model at an altitude h=15,000 feet and Mach numbers M=0.6 and 0.9. The 
analytical model for C nQ is constructed at p = 0°, 20°; stabilator deflections 8h = 10.5° 

and -24°; 8a=+ 25°; and 8r = + 30°. The analytical models are functions of a from 0° to 
90°; they are defined in Tables 8.1 a,b,c,d and 8.2 a,b,c,d. The analytical formulae are 
presented in Tables 8.3 a,b,c,d. Comparisons of the analytical models with the 
corresponding wind-tunnel model data are shown in Figures 8.1 to 8.7. The roll and 
yaw rate derivatives C„ p and C„ r and the sideslip derivative C n p are given in Figure 8.7. 

The analytical models are also given in the computer code listing contained in Appendix 
C. 


Table 8.1a Definitions of C n Analytical Models at M=0.6,8h=10.5° 


C no (a, J3=20°, M=0.6, 8a= 25°, 8r=-30”, h=15,000 ft ) 
C no (a, P=20°, M=0.6, 5a= 25°, 8r= 30°, h=15,000 ft ) 
C no (a, P=20°, M=0.6, 8a=-25”, 8r=-30”, h=15,000 ft ) 
C no (a, P=20°, M=0.6, 8a=-25”, Sr= 30°, h=15,000 ft ) 
C no (a, P=0° , M=0.6, 5a= 25°, 6r=-30°, h=15,000 ft ) 
C no (a, P=0° , M=0.6, 8a= 25°, 8r= 30”, h=15,000 ft ) 
C np (a, M=0.6, h=15,000 ft ) 

C nr (a, M=0.6, h=15,000 ft ) 

C np (M=0.6, h= 15,000 ft ) 


= C n XNXB2(a) 
= C n XXXB2(a) 
= C n NNXB2(a) 
= C n NXXB2(a) 
= CnXNXBO(a) 
= C n XXXB0(a) 
= C n P(a) 

= C n R(a) 

= C n B(a) 


Table 8.1b Definitions of C n Analytical Models at M=0.6,8h=-24° 


Cn 0 (oc, p=20°, M=0.6, 5a= 25”, 8r=-30”, h=15,000 ft ) 

C no (a, P=20°, M=0.6, 8a= 25°, 8r= 30°, h=15,000 ft ) 

C no (a, p=20°, M=0.6, 8a=-25”, 8r=-30”, h=15,000 ft ) 
C no (a, P=20°, M=0.6, Sa=-25”, 8r= 30°, h=15,000 ft ) 
C nfl (a, p=0° , M=0.6, 8a= 25°, 8r=-30«, h=15,000 ft ) 

C no (a, P=0« , M=0.6, 8a= 25”, 8r= 30°, h=15,000 ft ) 


= C n XNNB2(a) 
= C n XXNB2(a) 
= C n NNNB2(a) 
= C n NXNB2(a) 
= CnXNNBO(a) 
= C n XXNB0(a) 


Table 8.1c Definitions of Cn Analytical Models at M=0.9,8h=10.5° 


C no (oi, P=20°, M=0.9, 8a= 25”, 8r=-30”, h=15,000 ft ) 

C no (a, P=20°, M=0.9, 8a= 25°, 8r= 30°, h=15,000 ft ) 

C no (a, P=20°, M=0.9, 8a=-25", 8r=-30”, h=15,000 ft ) 
Cn 0 («»P= 20 °» M=0.9, 8a=-25°, 8r= 30”, h=15,000 ft ) 
C no («, P=0° , M=0.9, 8a= 25°, 8r=-30°, h=15,000 ft ) 

C no (a, P=0° , M=0.9, 8a= 25”, 8r= 30”, h=15,000 ft ) 


= C n XNX2(a) 
= C n XXX2(a) 
= C n NNX2(a) 
= C n NXX2(a) 
= CnXNXO(a) 
= C n XXX0(a) 
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Table 8.1d Definitions of Cn Analytical Models at M=0.9,8h=-24° 


C no (a»P=20°, M=0.9, 8a= 25«>, 8r=-30«, h=15,000 ft ) 
C no (a, P=20°, M=0.9, 8a= 25°, 8r= 30 H , h=15,000 ft ) 
C no (a, P=20°, M=0.9, 8a=-25°, 8r=-30°, h=15,000 ft ) 

C no (a, p=20°, M=0.9, 8a=-25°, 8r= 30°, h=I5,000 ft ) 

Cn 0 (a, P=0° , M=0.9, 8a= 25°, 8r=-30°, h=15,000 ft ) 

Cn 0 (a, 0=0° , M=0.9, 8a= 25°, 8r= 30°, h=15,000 ft ) 


= C n XNN2(a) 
= C„XXN2(a) 
= C n NNN2(a) 
= C n NXN2(a) 
= CnXNNO(a) 
= CnXXNO(a) 
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Table 8.2a Yawing Moment Coefficient Analytical Models: M = 0.6, Sh=10.5 


Sideslip 

Rudder 

Aileron 

Sr = -30 

8r = 30 

P = 20 

5a = 25 

CnXNXB2(a) 

CnXXXB2(a) 

5a = -25 

CnNNXB2(a) 

CnNXXB2(a) 

1 

6a = 25 

CnXNXBO(a) 

CnXXXBO(a) 

6a = -25 

-CnXXXBO(a) 

-CnXNXBO(a) 




Table 8.2c Yawing Moment Coefficient Analytical Models: M = 0.9, 8h=10.S 


Sideslip 

Rudder 
Aileron s. 

Sr = -30 

8r = 30 

P = 20 

6a = 25 

CnXNX2(ct) 

CnXXX2(a) 

8a = -25 

CnNNX2(ot) 

CnNXX2(a) 

p = o 

5a = 25 

CnXNXO(a) 

CnXXXO(a) 

8a = -25 

-CnXXXO(a) 

-CnXNXO(a) 


Table 8.2b Yawing Moment Coefficient Analytical Models: M = 0.6, 6h=-24 


Sideslip 

N. Rudder 

Aileron s. 

6r = -30 

6r = 30 

P = 20 

5a = 25 

CnXNNB2(cc) 

CnXXNB2(a) 

5a = -25 

CnNNNB2(a) 

CnNXNB2(a) 

xs 

II 

o 

6a = 25 

CnXN.NBO(a) 

CnXXNBO(a) 

5a = -25 

-CnXXNBO(a) 

-CnXNNBO(o) 


Table 8Jd Yawing Moment Coefficient Analytical Models: M = 0.9, 6h=-24 


Sideslip 

Rudder 

Aileron 

Sr = -30 

5r = 30 

P = 20 

8a = 25 

CnXNN2(a) 

CnXXN2(a) 

6a = -25 

CnNNN2(a) 

CnNXN2(a) 

p = o 

1 

6a = 25 

CnXNNO(a) 

CnXXNO(a) 

6a = -25 

-CnXXNO(a) 

-CnXNNO(a) 










































Table 8.3a Formulas for Cn Model at M=0.6, 8h=10.5° 


CnXNXB 2 (a°)=(. 06482/2. 75)tan* 1 (-(a°-20)3/14) 
+(.04937/2.75)tair 1 (-(a°-4 1)5/14) 
+(.053/2.75)tan- 1 ((a°-60)l/6) 

+ (. 005/2. 75)tarr 1 (-(a°-8)4) + .0211 

CnXXXB2(a°)=(. 077832/2. 75)tair 1 (-(oc°-27. 5)4/45) 

+(. 0744/2. 75)tair I ((a°-58. 5)4/25) 

+(. 02885/2. 75)tair 1 (-(a°-79)7/20) 
+(.006/2.75)tan- , (-(a°-43)) - .022 

CnNNXB 2 (a°)=(. 11977/2. 75)tan* 1 (-(oc°-28. 5)3.5/51) 
+(.04303/2. 75) tan’ 1 ((«°-58. 5)4/13) 

+ (. 02532/2. 75)tan _1 (-(ct 0 -74)2/5) + .01184 

CnNXXB2(a°)=(. 06132/2. 75)tan 1 (-(«°-31)I/10) 

+(. 05521/2. 75)tair , ((a°-55)2/9) 

+ (. 04659/2. 75) tarr 1 (-(a°-77. 5)4/15) - .03235 

CnXNX BO (a°)=(.02026/2.75)tair , (-(a°- 19)3/14) 

+(. 022/2. 75)tair 1 (-(a°-49)2/7) 
+(.03393/2.75)tarr 1 ((a°-73)l/7) 

+(. 002/2. 75)tan* 1 (-(a°- 14)2) 

+ (. 003/2. 75)tair 1 ((a°-76)2) + .02769 

CnX X XB 0(a°)=(. 00952/2. 75)tan- 1 ((cx°- 14)3/8) 

+(. 01056/2. 75)tair 1 ((a°-47)2/3) 
+(.01395/2.75)tan- l ((a°-67)9/14) 

+ (. 00899/2. 75) tan‘ 1 (-(a°-8 1)3/8)-. 01862 

CnNNXB0(a°)=-CnXXXB0(a°) 

CnNXXB0(a°)=-CnXNXB0(a°) 

CnP(a°) = (.075/Tt)tan-i((a°-17)5Tc/18) 
+(.04/7C)tan- 1 ((a°-50)10jc/18) 
+(.2/ic)tan* , (-(oc°-57)1007c/18) 
+(.13/7t)tan* 1 ((a°-62)1007t/18) 
+(.09/jt)tair 1 (-(a°-73)100Jc/18) 
+(.l/jt)tan‘ 1 ((a°-77)100jc/18)-.028 

CnR(a°) = (.16/7t)tan-»(-(a°-22)107c/18) 
+(.34/jr)tan’ 1 ((ct°-57)107c/18) 
+(-.l)exp((a°-78)l/10)-.09 

CnB(a°) = (.000001)[(12.7/Tc)tair 1 (-(a°-13) 10071/18) - 11.7] 


48 



Table 8.3b Formulas for C n Model at M=0.6, 8h=-24° 


CnXNNB2(cc°)=(.l 1857/2. 75)tan* 1 (-(a°-26. 5)4/45) 
+(.07065/2. 75)tair 1 ((a°-60. 5)4/25) 
+(.02518/2.75)tan 1 (-(a°-80)3/10) 
+(.005/2.75) tan' 1 (-(a°-18)2) + .020265 

CnXXNB2(a°)=(. 080916/2. 75)tan* 1 (-(a°-30)3/42) 

+(.056/2. 75)tan _1 ((ot°-62) 1/4) 

+(. 02085/2. 75)tan’ 1 (-(a°-79) 1/5) 
+(-.005/2.75)tan-i(-(a°-82)2) - .02221 

C nNNNB2(a°)=(. 120543/2. 75)tair 1 (-(a°-28)3/40) 

+(. 05707/2. 75)tan' 1 ((a°-55. 5)4/25) 

+ (. 03650/2. 75)tan* 1 (-(a°-78. 5)4/15) + .01202 

CnNXNB2(a°)=(.063768/2.75)tan’ 1 (-(a°-32)l/10) 

+(.04788/2. 75) tan* 1 ((a°-57. 5)6/25) 

+(. 03829/2. 75)tair 1 (-(a°-77. 5)4/15) - .03288 

Cn XNN B0(a°)=(. 02037/2. 75)tarr 1 (-(a°- 17.5)8/43) 
+(.00389/2.75)tan- , (-(a°-54.5)6/17) 

+(. 01 623/2. 75)tan- 1 ((cc°-69. 5)5/23) 

+ (. 002/2. 75)tan _1 (-(a°-13)2) + .02711 

CnXXNB0(a o )=(.00953/2.75)tan* 1 ((a o - 15)3/8) 

+(.00411/2.75)tan- 1 (-(a°-46. 5)8/15) 

+ (. 02222/2. 75)tan* 1 ((a°-71.5)4/25) - .01781 

CnNNNB0(a°)=-CnXXNB0(a°) 

CnNXNB0(a°)=-CnXNNB0(a°) 
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Tabic 8.3c Formulas for Cn Model at M=0.9, 8h=10.5 M 


CnXNX2(a°)=(.05428/2.75)tan' 1 (-(a°-16)3/14) 

+(. 05037/2. 75)tan- 1 (-(a°-36)5/14) 
+(-.003/2.75)tan' 1 ((a°-45)l/6) 

+ (.060/2.75)tan- 1 ((a°-50)l/6) 

+ (.005/2.75)tan- 1 (-(a°-8)4) + .0211 

CnXXX2(a°)=(.067832/2.75)tan' 1 (-(oc°-30)4/45) 
+(. 065/2. 75)tan' 1 (-(a°- 14)4/45) 

+ (-. 04/2. 75)tan* 1 (-(a°- 14)4/45) 

+ (. 0844/2. 75)tan' 1 ((a°-46)4/25) 

+ (.04085/2. 75)tan _1 (-(a 0 - 100)7/20) 

+ (. 006/2. 75) tan 1 ('(«°-35)) - .0352 

CnNNX2(a°)=(.05428/2.75)tan- 1 (-(a°-22)3/14) 

+ (.02037/2.75)tan' 1 (-(a°-36)5/14) 

+ (-.003/2. 75)tan- J ((a 0 -45) 1/6) 

+(. 042/2. 75) tan -1 ((a°-47) 1/6) 

+ (.020/2. 75) tan' 1 (-(a°-ll)4) + .0181 

CnNXX2(a°)=(.067832/2.75)tan- 1 (-(a°-35)4/45) 

+(. 070/2. 75)tan' 1 (-(a°- 14)4/45) 

+ (-.04/2.75)tan- 1 (-(a°-4)4/45) 

+ (. 0844/2. 75)tan' 1 ((a°-46)4/25) 

+ (. 04085/2. 75)tan' 1 (-(a°- 100)7/20) 

+ (. 006/2. 75)tan- 1 (-(a°-35)) - .0392 

CnXNX0(a°)=(.01026/2.75)tan' 1 (-(a°- 13)3/14) 

+ (-. 009/2. 75)tan' 1 (-(a°-40)2/7) 

+ (.010/2. 75)tan- 1 (-(a°-18)2/7) 

+ (-. 002/2. 75)tan' 1 (-(a°-30)2/7) 

+ (. 022/2. 75)tan' 1 (-(a°-49)2/7) 

+ (.02543/2. 75)tan‘ 1 ((a°-83) 1/7) 

+ (.002/2. 75)tan- 1 (-(a°-7)2) 

+ (. 003/2. 75)tan- 1 ((a°-76)2) + .02769 

CnXXX0(a°)=[(.01452/2.75)tan' 1 ((a°- 11)3/8) 

+ (-. 005/2. 75)tan' 1 ((a°-22)3/8) 

+ (. 01 156/2. 75)tan' 1 ((a°-41)2/3) 

+ (. 01205/2. 75)tan* 1 ((a°-67)9/14) 

+(. 00769/2. 75)tan' 1 (-(a°-81)3/8) 

- . 018621 ( 1 / 1 . 21 ) 

CnNNX0(a°)=-CnXXX0(a°) 

CnNXX0(a°)=-CnXNX0(a°) 
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Table 8.3d Formulas for Cn Model at M=0.9, 5h=-24° 

CnXNN2(a°)=(.05428/2.75)tair 1 (-(a°-15)3/14) 

+(. 05037/2. 75)tair 1 (-(a°-33)5/14) 
+(-.003/2.75)tair 1 ((a°-45)l/6) 
+(.060/2.75)tan* , ((a°-46)l/6) 
+(.005/2.75)tan-H-(a°-8)4)+.0241 

CnXXN2(a°)=(.11091/2.75)tair 1 (-(a°-25)3/42) 

+ (-. 025/2. 75)tan* 1 (-(a°-8)3/42) 

+(.05600/2.75)tan- 1 ((a°-48)l/4) 

+(.03385/2.75)tan- 1 (-(a°-100)l/5) 

+(-. 005/2. 75)tair 1 (-(a°-82)2) - .03125 

CnNNN2(a°)=CnXNN2(a°) 

CnNXN2(a°)=(. 067832/2. 75)tair 1 (-(a°-32)4/45) 

+ (.063/2.75)tan-i(-(a°-14)4/45) 

+(-.04/2.75)tan-i(-(a°-4)4/45) 

+(.0844/2.75)tan* 1 ((a°-46)4/25) 

+(. 04085/2. 75)tair 1 (-(a°-90)7/20) 

+(. 006/2. 75)tair 1 (-(a°-35)) 

-.0352 


CnX NN0(a°)=(. 01637/2. 75)tan* , (-(a°- 10)8/43) 

+ (.00559/2.75)tan- 1 (-(a°- 180)6/17) 

+ (. 01623/2. 75)tair 1 (-(a°- 100)5/23) 
+(.002/2.75)tan* 1 (-(oc°-13)2) 

+.0271 

CnXXN0(a°)=(.01253/2.75)tan- 1 ((a°-12)3/8) 

+(-. 002/2. 75)tarr 1 ((a°-22)3/8) 

+ (.00411/2.75)tan- 1 (-(a°-46.5)8/15) 

+ (. 02222/2. 75)tan* 1 ((«°-81. 5)4/25) - .01161 

CnNNN0(a°)=-CnXXN0(a°) 

CnNXN0(a°)=-CnXNN0(a°) 
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Figure 8.1: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient C, 
for h=15,000 feet , M=0.6, 8h= 10.5° and p=2(R 
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Figure 8.2: Comparison of Wind-Tunnel and Analytical Yawing 
for h= 15,000 feet , M=0.9, 8h= 10.5° and p=20°. 
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Figure 8.3: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient C^ 
for h= 15,000 feet , M=0.6, 8h= -24° and p=20°. 
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Figure 8.4: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient C, 
for h= 15,000 feet , M=0.9, 8h= -24° and p=20°. 
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Figure 8.5: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient C. 


"0 


for h= 15,000 feet , M=0.6, 8a= 25° and (3=0°. 


56 




CNO CNO 


Analytical Model of 


Wind-Tunnel 


CNXNXO(tx) 



8h= 10.5°, Sr= -30° 


CNXNNO(a) 



Analytical 


CNXXXO(a) 



6h= 10.5°, 8r= 30° 


CNXXNO(a) 

o.oi — 



-0.03- 

*°‘ M ■ I ~~T 1 — I V ' | T 

-«.0 -20.0 0.0 20.0 «. 0 60.0 00.0 100.0 120.0 
RLPDE6 


8h= -240, 5r= .300 


8h= -24°, 5r= 30° 


Figure 8.6: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient 
for h= 15,000 feet , M=0.9, 8a= 25° and p=0°. 
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M=0.6, h=15,000 feet 
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CnB(a) [ C np (a)] 



M=0.6, h=15,000 feet 

Figure 8.7: Comparison of Wind-Tunnel and Analytical Yawing Moment Coefficient 
Derivatives C,^ , C,^ and C 0 p for h=15,000 feet and M=0.6. 
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9. Time History Comparison of u,v,w,p,4 and f : Mach = 0.6 


The simulation wind-tunnel model of [1] was flown in NASA's simulator by a pilot 
to generate some basic maneuvers at 0.6 Mach numbers such as pitch-ups, 360° loaded and 
unloaded rolls, turn reversals, split S's and level turns. That simulator data is used here to 
check the validity of the 6 DOF analytical model. The accelerations 

u, w, q, v, p, r 

are computed for the analytical model using the states and controls from the piloted 
simulated maneuvers. Comparisons with the accelerations from the wind-tunnel data 
model are shown below in Figures 9. 1-9.6. 
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U.DOT 


Analytical Model Simulation: Pitch Up Maneuver @ M=0.6 (Run 1, 6 October 1987) 
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Figure 9.1: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Pitch Up Maneuver @ M=0.6 (Run 1, 6 October 1987). 
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Analytical Model Simulation: 360° Loaded Roll (Trim Power) Maneuver @ M=0.6 (Run 3, 6 October 1987) 
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Figure 9.2: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
360° Loaded Roll (Trim Power) Maneuver @ M=0.6 (Run 3, 6 October 1987). 


u, w, q, v, p, r 




V.DOT . U.DOT 


9 (deg/sec 2 ) 


Analytical Model Simulation: Turn Reversal Maneuver @ M=0.6 (Run 4, 6 October 1987) 
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Figure 9.3: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Turn Reversal Maneuver @ M=0.6 (Run 4, 6 October 1987). 
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Analytical Model Simulation: 360° Loaded Roll (AB) Maneuver @ M=0.6 (Run 5, 6 October 1987) 
u (ft/sec 2 ) w (ft/sec 2 ) Q (deg/sec 2 ) 



Wind-Tunnel Analytical 

Figure 9.4: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: u » q, v, p, r 

NASA's 360° Loaded Roll (AB) Maneuver @ M=0.6 (Run 5, 6 October 1987). 
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Figure 9.5: Comparison of Derivatives Generated by W; 

Split S Maneuver @ M=0.6 (Run 6, 6 Octot 
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Analytical Model Simulation: Level Turn Maneuver @ M=0.6 (Run 7, 6 October 1987) 
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Figure 9.6: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Level Turn Maneuver @ M=0.6 (Run 7, 6 October 1987). 
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10. Time History Comparison of ii,v,w,p,q and i : Mach = 0.9 


The simulation wind-tunnel model of [1] was flown in NASA's simulator by a pilot 
to generate some basic maneuvers at 0.9 Mach numbers such as pitch-ups, 360° loaded and 
unloaded rolls, turn reversals, split S's and level turns. That simulator data is used here to 
check the validity of the 6 DOF analytical model. The accelerations 

u, w, q, v, p, r 

are computed for the analytical model using the states and controls from the piloted 
simulated maneuvers. Comparisons with the accelerations from the wind-tunnel data 
model are shown below in Figures 10.1-10.7. 

The piloted simulated maneuvers comparison herein shows that the angular pitch 
accelerations from the wind-tunnel data and the analytical model have about the same shape 
but at times have a fairly large distance between them. We show in Appendix E that this is 
due to a small error in fit being multiplied by a large dynamic pressure at Mach 0.9. 
Therein we present some details from Run 5 (Figure 10.3) which is a turn reversal 
maneuver to show that the differences are due to a small difference of about 0.006 or less 
in the values of C mQ (t). As can be seen from the modeling fits shown in Figures 5. 1-5.6 

modeling errors of this magnitude are present in C mQ at all Mach numbers. We found that 

the largest differences are equivalent to approximately a half degree change in stabilator 
deflection angle. 
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Analytical Model Simulation: Pitch Up Maneuver @ M=0.9 (Run 3, 10 October 1987) 
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Figure 10 . 1 : Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model- w, q, V, p, r 

Pitch Up Maneuver @ M=0.9 (Run 3, 10 October 1987). 
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Analytical Model Simulation: Turn Reversal Maneuver @ M=0.9 (Run 5, 10 October 1987) 
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Figure 10.3: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model* W, q, v, p, r 

Turn Reversal Maneuver @ M=0.9 (Run 5, 10 October 1987). 







Analytical Model Simulation: 360° Unloaded Roll (AB) Maneuver @ M=0.9 (Run 9, 10 October 1987) 
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Figure 10.4: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
360° Unloaded Roll (AB) Maneuver @ M=0.9 (Run 9, 10 October 1987). 
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Analytical Model Simulation: Split S Maneuver @ M=0.9 (Run 6, 10 October 1987) 
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Figure 10.6: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Level Turn Maneuver @ M=0.9 (Run 7, 10 October 1987). 
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Analytical Model Simulation: 360° Unloaded Roll (MIL PWR) Maneuver @ M=0.9 (Run 8, 10 October 1987) 
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Figure 10.7: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
360° Unloaded Roll (MIL PWR) Maneuver @ M=0.9 (Run 8, 10 October 1987). 
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11. Time History Comparison of u,v,w,p,q and r : Mach = 0.3 


The simulation wind-tunnel model of [ 1] was flown in NASA’s simulator by a pilot 
to generate some basic maneuvers at 0.3 Mach numbers such as pitch-ups, 360° loaded and 
unloaded rolls, turn reversals, split S's and level turns. That simulator data is used here to 
check the validity of the 6 DOF analytical model. The accelerations 

li, w, q, v, p, r 

are computed for the analytical model using the states and controls from the piloted 
simulated maneuvers. Comparisons with the accelerations from the wind-tunnel data 
model are shown below in Figures 1 1.1-1 1.7. 
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Analytical Model Simulation: Pitch Up Maneuver @ M=0.3 (Run 11, 6 October 1987) 
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Figure 11.1: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Pitch Up Maneuver @ M=0.3 (Run 1 1, 6 October 1987). 
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Analytical Model Simulation: 360° Loaded Roll Maneuver @ M=0.3 (Run 13, 6 October 1987) 
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Figure 11.2: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
360° Loaded Roll Maneuver @ M=0.3 (Run 13, 6 October 1987). 
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Figure 11.3: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Turn Reversal Maneuver® M=0.3 (Run 15, 6 October 1987). 
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Analytical Model Simulation: 360° Unloaded Roll Maneuver @ M=0.3 (Run 21, 6 October 1987) 
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Figure 1 1.4: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
360° Unloaded Roll Maneuver @ M=0.3 (Run 21, 6 October 1987). 
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Analytical Model Simulation: Split S Maneuver @ M=0.3 (Run 16, 6 October 1987) 
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Figure 11.5: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Split S Maneuver @ M=0.3 (Run 16, 6 October 1987). 
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Analytical Model Simulation: Level Turn (MAX) Maneuver @ M=0.3 (Run 20, 6 October 1987) 
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Figure 11.6: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model- U, W, q, v, p, r 
Level Turn (MAX) Maneuver @ M=0.3 (Run 20, 6 October 1987). 
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Analytical Model Simulation: Split S Maneuver @ M=0.3 (Run 18, 6 October 1987) 
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Figure 11.7: Comparison of Derivatives Generated by Wind-Tunnel and Analytical Model: 
Split S Maneuver @ M=0.3 (Run 18, 6 October 1987). 




12. SUMMARY 


A six degrees of freedom (6 DOF) analytical aerodynamic model is derived from a 
high angle-of-attack combat airplane wind-tunnel model, [1]. The derivation considered 
the altitude-Mach flight envelope centered at an altitude h = 15,000 feet and a Mach number 
M = 0.6. Wind-tunnel data ranging from 0.3 Mach to 0.9 Mach was used in developing 
the analytical models. The derived analytical models of the aerodynamic derivatives are 
nonlinear functions of alpha with all other states and control variables fixed. The nonlinear 
functions are parameterized with respect to sideslip, Mach number, roll, pitch and yaw 
rates and aileron deflection, rudder deflection and stabilator deflection. The lift and 
pitching moment coefficients have unsteady flow parts due to the time rate of change of 
angle of attack (alpha dot). The effects of leading edge flap, trailing edge flap, speed 
brake, landing gear, etc was not consided. Interpolation is required between the 
parameterized nonlinear functions. 


Formulae for the analytical models are shown to compare well with their fits of the 
wind-tunnel data. The piloted simulated maneuvers comparison of Chapters 9-11 in which 
the analytical model is compared with the wind-tunnel data show that (1) the analytical 
model is a good representation of the wind-tunnel at Mach 0.6, (2) the longitudinal part of 
the analytical model is good for the Mach number range 0.3 to 0.9 and (3) the lateral part is 
good for Mach numbers between 0.6 to 0.9. Analytical models of the rolling moment 
coefficient were not derived using 0.3 Mach wind-tunnel data. The piloted simulated 
maneuvers comparison of Chapter 1 1 indicates that analytical models of the rolling moment 
coefficient should be fit at Mach 0.3 in order to better represent the wind-tunnel model 
there. The piloted simulated maneuvers comparison of Chapter 10 shows that the angular 
pitch accelerations from the wind-tunnel data and the analytical model have about the same 
shape but at times have a fairly large distance between them; this is due to a small C mQ error 

of about 0.006 in fit being multiplied by a large dynamic pressure at Mach 0.9; see 
appendix E. We found that the largest differences are equivalent to approximately a half 
degree change in stabilator deflection. Such a difference should not be critical in analysis 
studies. 


The results in this report indicate that the analytical model is a good representation 
of the wind-tunnel model for flight analysis in the altitude-Mach flight envelope centered at 
an altitude h = 15,000 feet and a Mach number M = 0.6. The storage requirement of the 
analytical model is about one tenth that of the wind-tunnel model and it runs twice as fast. 
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APPENDIX A 

Computer Code for Equations of Motion 
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== 

COMAN 

THIS PROGRAM IS USED TO COMPUTE U. DOT, V. DOT 

= = 

== 


W. DOT, P. DOT 

, Q. DOT, R. DOT FOR THE HARV ANALY- 

== 

== 


TICAL MODEL 


== 

=s 

INPUTS: 

FLIGHT DATA 

FILE 

ass 

*=« 


T 

TIME 

== 

== 


MACH 

MACH NUMBER 

== 

== 


HAB 

ALTITUDE 

== 

ss 


QBAR 

DYNAMIC PRESSURE 

== 

== 


ALPDEG 

ANGEL OF ATTACK 

== 

== 


BETADEG 

SIDESLIP ANGEL 

== 

== 


PH ID 

EULER BANK ANGEL 

== 

== 


THETAD 

EULER PITCH ANGEL 

== 

== 


PS ID 

EULER YAW ANGEL 

== 

== 


P 

AIRCRAFT X-BODY AXIS ROLL RATE 

== 

— 


Q 

AIRCRAFT Y-BODY AXIS PITCH RATE 

== 



R 

AIRCRAFT Z-BODY AXIS YAW RATE 

== 

== 


DH 

STABILATOR DEFLECTION 

== 

== 


DA 

AILERON DEFLECTION 

ss 



DR 

RUDDER DEFLECTION 

ss 

== 

OUTPUTS: 

UDT 

TIME DERIVATIVE OF U 

sss 

== 


VDT 

TIME DERIVATIVE OF V 

== 

== 


WDT 

TIME DERIVATIVE OF W 

== 

== 


PDT 

TIME DERIVATIVE OF P 

ss 

== 


QDT 

TIME DERIVATIVE OF Q 

== 

== 


RDT 

TIME DERIVATIVE OF R 

== 


== AUTHER J I CHANG CAO 

== GRADUATE RESEARCH ASSISTANT 

== SCHOOL OF AEROSPACE ENGINEERING 

== GEORGIA INSTITUTE OF TECHNOLOGY 

== ATLANTA, GEORG I A 30332 

== DATA JUNE, 1990 


PROGRAM COMLAT 
PARAMETER (N=129) 

IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

REAL *8 MACH .MASS ,NZ , I AS ,HAB 
REAL *8 AU(N) ,AV(N) ,AW(N) ,AP(N) ,AQ(N) ,AR(N) 
real *8 APHID(N) ,ATHETA(N) ,APSID(N) ,AALP(N) ,ABET(N) 
real *8 TIME(N) ,AMACH(N) ,ahab(n) ,aqbar(n) ,avtotal(n) 
real *8 ADH(N) ,ADA(N) ,ADR(N) 

REAL *8 UDT(N) ,VDT(N) ,WDT(N) ,PDT(N) ,QDT(N) ,RDT(N) 
real *8 ATNETL (N) .ATNETR(N) ,lxe ,lye ,lze 
character* 8 o header 

CHARACTER AA1*15, AA2*15, AA3* 1 5 , AA4* 1 5 , AA5*15,BB*50 
OPEN (UNIT=5,FILE='MACH03') 

OPEN (UNIT= 6 ,FILE=’DATA_ANAL') 


GENERAL CONSTANTS 


G =32.174 

DTR =ACOS (- 1 .) / 1 80 . 


AIRCRAFT CONSTANTS 


MASS =1035-308 

XI = 23000 . 

Y I =151293. 

zi = 169945 . 
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XZI =-2971. 

B =37.42 

CBAR =11.52 

S =400. 

XL =-3.56/12. 

VL =0. 

ZL =2.8/12. 

THZ =0.0 

LXE =-232.5/12. 

LYE =0.0 

LZE =2.8/12. 

CSTAR =X I *Z I / (X I *1 1 -XZ I **2) 

C4 1 =CSTAR*XZ I * (Z l+X I -Y I ) / (X I *Z I ) 

C42 =CSTAR* (Z I * (Y I -Z I ) -XZ I **2) / (X I *Z I ) 
C43 =CSTAR*XZ I /X I 

C5 1 = (Z t -X I ) /Y I 

C52 =XZ 1 /Y I 

C6 1 =CSTAR* (X I * (X I -Y I ) +XZ I **2) / (X I *Z I ) 

C62 =CSTAR*XZ I * (Y I -Z I -X I ) / (X I *Z I ) 

C 63 =CSTAR*XZ I /Z I 


FLIGHT CONDITION 


READ (5,5) NUMBER 
DO 99 1=2,13 

READ (5,15) AA1 , AA2 ,AA3 ,AA4 ,AA5 
99 CONTINUE 

READ (5,25) BB 
5 FORMAT (14) 

15 FORMAT (5A 1 5) 

25 FORMAT (A50) 


DO 400 1=1, N 

READ FLIGHT DATA FILE 

READ (5,10) T .MACH , HAB , QBAR .ALPDEG .BETADEG ,PSID .THETAD 
READ (5,10) PHID , PDEG , QDEG , RDEG ,VTOTAL,AX , AY ,GZ 
READ (5,10) DLADEG,DLADDEG,DLSDEG,DLHDDEG,DLNDEG,DLNDDEG, 

1 DLFDEG . DLFDDEG 

READ (5,20) DLRDEG , DLRDDEG .TNETL .TNETR 


ft ft ft ft ft ft ft ft ft ft ft ft ft * ft ft ft ft ft ft ft ft ft * * * ft * ft ft ft ft j't :'t )'t A :'c )'c ;V )'t :'c * :V ,'c i; i: * ;'c :'c ;'c ;'c )V ;V ;'c )'t ;V 


10 

20 


T 

— 

TIME 

MACH 

— 

MACH NUMBER 

HAB 

— 

ALTITUDE 

QBAR 

— 

DYNAMIC PRESSURE 

ALPDEG 

— 

ANGEL OF ATTACK 

BETADEG 

— 

SIDESLIP ANGEL 

PS ID 

— 

EULER BANK ANGEL 

THETAD 

— 

EULER PITCH ANGEL 

PHID 

— 

EULER YAW ANGEL 

PDEG 

— 

AIRCRAFT X-BODY AXIS ROLL RATE 

QDEG 

— 

AIRCRAFT Y-BODY AXIS PITCH RATE 

RDEG 

— 

AIRCRAFT Z-BODY AXIS YAW RATE 

VTOTAL 

— 

VELOCITY 

DLADEG 

— 

AILERON DEFLECTION (AVERAGE) 

DLRDEG 

— 

RUDDER DEFLECTION (AVERAGE) 

DLSDEG 

— 

STABILATOR DEFLECTION (AVERAGE) 

TNETL 

— 

THRUST OF LEFT ENGINE 

TNETR 

— 

THRUST OF RIGHT ENGINE 


it * #c it sV )'c i'c ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft 

FORMAT 

(8E10 

.4) 

FORMAT 

(4E10 

.4) 


i’c jV )'c A sV sV :'c A :'c :V ;V i't 5'; 


= VTOTAL*DCOS (ALPDEG*DTR) *DCOS (BETADEG*DTR) 
= VTOTAL*DSIN (BETADEG*DTR) 


U 

V 
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w = vtotal*dsin (ALPDEG*DTR) *dcos (betadeg*dtr) 


AVTOTAL (1) 

i = VTOTAL 

AU ( 1 ) = U 


AV ( 1 ) = V 


AW (1 ) = W 


AP ( 1 ) = PDEG 

AQ ( 1 ) = QDEG 

AR ( 1 ) = RDEG 

APHID (1) 

=PH 1 D 

ATHETA(I) 

=THETAD 

APS 1 D ( 1 ) 

=PS 1 D 

AALP(I) 

=ALPDEG 

ABET { 1 ) 

=BETADEG 

AHAB(I) 

=HAB 

TIME (1) 

=T 

AQBAR(I) 

=QBAR 

AMACH (1) 

=MACH 

ADH ( 1 ) 

=DLSDEG 

ADA ( 1 ) 

=DLADEG 

ADR ( 1 ) 

=DLRDEG 

ATNETL (1) 

=TNETL 

ATNETR (1) 

=TNETR 


400 CONTINUE 

DO 410 1=1 ,N 

B02VT = B/(2*AVT0TAL(I)) 

C02VT = CBAR/(2*AVT0TAL (!)) 

ALPHA =AALP ( I ) 

BETA =ABET ( I ) 

MACH =AMACH ( I ) 

ALT =AHAB ( I ) 

DH =ADH ( I ) 

DA =ADA(l)/2. 

DR =ADR ( I ) 

QBAR =AQBAR ( I ) 

TNETL =ATNETL ( I ) 

TNETR =ATNETR ( I ) 

P =AP ( I ) *DTR 
Q =AQ (I) *DTR 
R =AR ( I ) *DTR 

CALL NAER0C1 (CLO ,CLQ .CLAD.CMO , CMQ , CMAD , 

* CDO , 

=F ( 

* ALPHA, MACH, DH, ALT) 

CALL NAER0C2 (CYO ,C10 ,CNO ,CYB ,C1B ,CNB , 

* CYR , C 1R , CNR , CYP ,C1P ,CNP , 

=F ( 

* ALPHA .BETA .MACH , DA , DR , DH ,P ,R .ALT ) 


U =AU ( I ) 

V =AV ( I ) 

W =AW ( I ) 


PH ID =APH I D ( I ) *DTR 
THET =ATHETA(I)*DTR 
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PS ID =APS I D ( I ) *DTR 


AA =ALPHA*DTR 

BT =BETA*DTR 

THX = (TNETL+TNETR) *COS (1 . 98*DTR) 

THY = TNETR*SIN(-1.98*DTR)+TNETL*SIN(1.98*DTR) 

FD = QBAR*S*CDO/MASS 

CB = QBAR*S*C02VT*CLAD/ (U*U+W*W) /MASS 

BQ = 1+CB*DC0S (AA) *U-CB*DCOS (AA) *W* (CB*DS I N (AA) *U) 

* / (1+CB*DS I N (AA) *W) 

FU = R*V-Q*W-G*DS I N (THET) -FD*DCOS (AA) +THX/MASS 

* +QBAR*S*DS 1 N ( AA ) * (CL0+C02VT*CLQ*Q) /MASS 

FW = Q*U-P*V+G*DCOS (THET) *DCOS (PHID) -FD*DSIN (AA) 

* -QBAR*S*DCOS (AA) * (CL0+C02VT*CLQ*Q) /MASS+THZ/MASS 
WDT (I) = (FW+CB*DCOS (AA) *W*FU/ (1+CB*DS I N (AA) *W) ) /BQ 

UDT ( I ) = (FU+CB*DS I N (AA) *U*WDT (I ) ) / ( 1+CB*DS I N (AA) *W) 

DALFA = (U*WDT ( I ) -W*UDT ( I ) ) / (U*U+W*W) 

CL = CLO+C02VT* (CLQ*Q+CLAD*DALFA) 

CM = CMO+C02VT* (CMQ*Q+CMAD*DALFA) 

Cl = C10+C1B*BT+B02VT*(C1P*P+C1R*R) 

CY = CYO+CYB*BT+B02VT* (CYP*P+CYR*R) 

CN = CN0+CNB*BT+B02VT*(CNP*P+CNR*R) 

FL = QBAR*S*CL/MASS 

FX = -FD*DCOS(AA)+FL*DSIN(AA) 

FY = QBAR*S*CY/MASS 

FZ = -FD*DS I N (AA) -FL*DCOS (AA) 

FP = QBAR*S*B*C1/XI+MASS*(YL*FZ-ZL*FY)/XI 

FQ = QBAR*S*CBAR*CM/YI+MASS*(ZL*FX-XL*FZ) /Y I 

FR = QBAR*S*B*CN/ZI+MASS*(XL*FY-YL*FX)/ZI 

VDT(I) = P*W-R*U+G*COS (THET) *S IN (PHID) +FY+THY/MASS 

YY1 = Citl*P*Q+CL2*Q*R+Ci43*FR+CSTAR*FP+Ci*3*LXE*THY/ZI 

1 -c43*LYE*THX/Z I 

1 +CSTAR* (LYE*THZ-LZE*THY) /X I 

YY2 = C51*P*R+C52* (R*R-P*P) +FQ+ (LZE*THX-LXE*THZ) /Y I 

YY3 = C6l*P*Q+C62*Q*R+C63*FP+CSTAR*FR+C63/X I * (LYE*THZ-LZE*THY) 

1 +cstar/zi*lxe*thy 

1 -CSTAR*LYE*THX/ZI 

PDT(I) =YY1/DTR 
QDT(I) =YY2/DTR 
RDT(I) =YY3/DTR 

LIO CONTINUE 

M=N-2 

DO 710 1=1, M 

WRITE (6,910) TIME ( I ) , UDT ( I ) ,VDT(I) , WDT ( I ) , PDT ( I ) ,QDT(I), 

1 RDT(I) 

710 CONTINUE 
910 FORMAT (7E10. 4) 

STOP 

END 
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APPENDIX B 

Computer Code for Longitudinal Analytical Model 
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non n ooono r> o ooooooooooooooooonooooooooooonooooo 


PURPOSE: THIS SUBROUTINE WILL BE CALLED TO CALCULATE 7 AERODY- 
NAMIC COEFFICIENTS WHICH ARE THE OUTPUTS OF THIS SUB- 
ROUT I NE . (LONG I TUD I NAL COEFFICIENTS) 


INPUTS: 

ALPDEG ANGLE OF ATTACK 

MACH MACH NUMBER 

ALT ALTITUDE 

DH STABILATOR DEFLECTION 

OUTPUTS: 

CLO , CLQ .CLAD 
CMO , CMQ , CMAD 
CDO 

THOSE COEFFICIENTS ARE USED IN THE FOLLOWING FORMULAS 

CL = CL0+C02VT*(CLQ*Q+CLAD*ALPDEG) 

CM = CM0+C02VT* (CMQ*Q+CMAD*ALPDEG) 

CD = CDO 

CL LIFT FORCE COEFFICIENT ALONG Z_WIND AXIS 

CM PITCH MOMENT COEFFICIENT ABOUT Y_WIND AXIS 

CD SIDE FORCE COEFFICIENT ABOUT X_WIND AXIS 

CREF REFERENCE WINGSPAN 

V AIRCRAFT TOTAL AIRSPEED 

AUTHER: J I CHANG CAO 

GRADUATE RESEARCH ASSISTANT 
GEORGIA INSTITUTE OF TECHNOLOGY 
ATLANTA, GA30332 


(DEG) 

(FEET) 

(DEG) 


SUBROUTINE NAEROC 


* 


(CLO , CLQ .CLAD, CMO , CMQ 
CDO , 

=F ( 

ALPDEG, MACH, DH, ALT) 


, CMAD, 


IMPLICIT REAL (C) 

REAL ALFA, MACH, DH, ALT. ALPDEG, A, 
DHN.DHX.PI 


CONVERSION OF ALFA 


IF (DH.LT.-2A.) DH=-24. 

IF (DH.GT.10.5) DH= 10.5 


A = ALPDEG 

PI =ACOS (- 1 . ) 

SI - 2.75 


EXTREMAL VALUES OF CONTROL SURFACES DEFLECTION 


DHN = -24 . 

DHX = 10.5 

c 

C COMPUTATION OF COEFFICIENTS 

c 

CL0X6 = 0.86/SI*ATAN (- (A+5.) jVI ./100.) 

* +2.19/SI*ATAN((A-5.)*l./7-) 

* +0.90/SI*ATAN ((A-24.) *1 ./I 7 -) 

* +1 .71/SI*ATAN (- (A-53-) *2./25.) 
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c 


c 


c 


c 


c 


c 


c 


c 


* +.4l/SI*ATAN (- (A-70.) *2 ,/7 .) 

* -.1+.05 

CL0N6 = 1 .06/SI*ATAN (- (A+5.) *1 ./100.) 

* + 1 . 79 /SI>vatan ( (A-5 •) *1 ./7.) 

* +2 . 50/S I AATAN ( (A- 1 5 .) *1 ./22 .) 

A +2 .71/S I AATAN (- (A -59 . ) *1 ./50.) 

* + 1 . 2 1 /S I *ATAN (- (A-70.) *1 ./20.) 

* -. 80+. 08 

CMOX 6 = 0.26/SI*ATAN (- (A-5.) *1 ./10.) 

* -0.39/SIAATAN ( (A-l .) *1 ./ 8 .) 

* +0.80/SI*ATAN((A-5.)*l./13.) 

* +0 . 70/S I AATAN (- (A- 10 .) *1 ./65.) 

* + 1 . 20/S I *AT AN ( (A-49 . ) * 1 . / 1 5 . ) 

* +2.1 0/S I AATAN (- (A-69 . ) * 1 / 1 5 . ) 

* -0.45/S I *ATAN(- (A-77.) *1 ./ 2 .) 

A -. 1 * 08+. 01 

CM006 = 0.36/SI*ATAN(-(A-5.)*l./30.) 

* -0 .29/SI AATAN ( (a - 1 . ) * 1 . / 1 5 . ) 

A +0 .90/S I *ATAN ( (A~5.) *1 ./35*) 

a + 0 . 80 /S 1 aatan (- (a-48.) ai V 75 .) 

* + 0 . 90 /S I AATAN ( (A - 52 . ) * 1 . / 1 0 . ) 

* + 2 . 10 /SI AATAN (- (A- 69 .) *1/15-) 

A -0.45/S I AATAN (- (A-77-) *1 ./2.) 

* -.1*5-- 007 

CM0Z6 = 0.26/SIAATAN(-(A-5.)A1./60.) 

* -0 . 39/S I AATAN ( (A - 1 .) A 1 . / 1 4 . 0) 

* +0 . 80/S I AATAN { (A-5 . ) A 1 . /42 . ) 

A +0.80/S I AATAN (- (A-20.) *1 ./55-) 

A + 1 . 80/S I AATAN ( ( A-65 • ) * 1 . /60 . ) 

* +2 . 40/S I AATAN (- (A- 69 • ) * 1 /20 . ) 

A -0 .45/SI AATAN (- (A-79 • ) * 1 . /2 . ) 

* -. 138-. 02 


CM0X26 = 0.26/SI *ATAN (- (A-2 .) *1 ./ 1 0 .) 

A -0.39/SIAATAN((A-1.)Al./]0.) 

a + 1 .00/SI AATAN ( (A-3 .) *1-/11.) 

* +0. 85/S I AATAN (- (A-7.) *1 ./20 .) 

* +1 .35/S I *ATAN ( (A-5 1 •) *1 -/1 9.) 

A + 2 . 20 /SI AATAN {- (A- 69 .) *1/15.) 

A -0.45/S I *ATAN (-(A-77.) A1./2.) 

* -. 3-. 01 

CMOX 56 = 0 . 26 /S I AATAN (- (A-5. ) * 1 ./ 10 .) 

* -0 . 39/S I AATAN ( (A- 1 . ) * 1 . / 1 2 . ) 

* + 0 . 90 /SIAATAN ((A-5.) * 1 ./15.) 

* + 0 . 85 /S I AATAN (- (A- 10 . ) * 1 ./30 •) 

a + 1 . 35 /s I aatan ((A-49.)*! . / 1 9 . ) 

A +2 . 20/S I *ATAN (- (A- 69 .) *1/15 •) 

A -0.45/S I AATAN (- (A-77.) A] */2 .) 

* -.348-. 02 

CM0N6 = 0.26/SI*ATAN(-(A-5.)*l./60.) 

A - 0 . 39 /S I *ATAN ( (A- 1 . ) A ] . / 30 .) 

* +0 . 80/S I *ATAN ( (A-5 . ) * 1 . /45 . ) 

A +0 . 80/S I AATAN (- (A- 1 0 . ) * 1 . /65 . ) 

ft +1 .80/S I AATAN ( (A-49.) *1 ./45.) 

A +2 . 80/S I AATAN (- (A-69-) *1/23-) 

A -0. 45/S I AATAN (-(A-79.) * 1./2.) 

* -. 138-. 01 


CMON 56 = 0 . 26 /S I *ATAN (- (A~5 ■ ) * 1 »/30 •) 
A -0.39/SIAATAN ((A-l .) *1 ./30.0) 
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c 


c 


c 


c 


c 


c 


* +1 . 20/S I *ATAN ( (A-5 - ) * 1 - AO.) 

* +0 . 60/S I SATAN (- (A- 8 . ) * 1 . /23 . ) 

* +1.30/SI*ATAN(-(A-60.)*l./65.) 

* +2 . 80/S I ftATAN ( (A-72 .) *1 ./55 •) 

* +2. 30/S I *ATAN (- (A-73-) *1/19-) 

* -0.45/S I *ATAN (- (A-77.) *1 ./ 2.) 

* -. 168-. 02 

CM 0 N 6 = 0 . 26 /S I SATAN (-(A- 5 JS 1 ./ 60 .) 
s -0. 39/S I *ATAN ( (A— 1 .) * 1 ./ 30 .) 

s + 0 . 80/S I SATAN ( (A-5 . ) * 1 . A5 • ) 

s +0.80/SISATAN (- (A- 10 .) si ./65 .) 

* + 1 .80/S I SATAN ( (A— 5 1 •) *1 . A5») 

* + 2 . 80 /S I SATAN (- (A- 69 .) * 1 / 23 .) 

* -0.45/S I *ATAN (- (A- 79 .) *1 ./ 2 .) 

s -. 138-. 01 

it it ic :'c * it : it i; it it * ;'c :'t it i< it :’c it * it it it ic ic it ic it jSt j'c j'c )'c j'c j'c j'c j'c * j'c it * j'c * j'c j'c * A A j’c ic j'c j'c j’c i< it j’c * j'c j'c 5V j’c j’c j'c j'c j'c 

CDOX = oAo/sisatan ((AS78./8o.+7.)si ./30.) 

s +o. 6 o/sisatan (- (AS78./80.+2.) si ./8.) 

* -0.30/S I SATAN (- (AS 78 ./ 8 O.+ 5 .) * 1 ./90.) 

s -0.20/S I SATAN (- (AS 78 ./ 80 .- 6 .) *1 ./5.) 

* +1 -95/S I SATAN ( (AS 78 ./ 8 O.- 28 .) *1 ./15.) 

* +2.20/S|SATAN((AS78./80.-58.)*l .AO.) 

* +1 .40/SI SATAN (- (AS 78 ./ 80 .- 73 .) SI ./ 30 .) 

s +2. 30/S I SATAN (- (AS 78 ./ 8 O.- 138.) *1/20.) 

* -.147 

CDOZ « (0.60/S I SATAN ((A*77./80.+6.) *1 ./30.) 

* +0. 60/SI SATAN (- (A*77./80. + l .) *1 ./ 8 .) 

* -0.30/S I SATAN (- (A*77 ./80 .+!».) *1 ./90.) 

* -0.20/S I SATAN (- (AS 77 ./80 . - 7 .) *1 ./TO.) 

* +1 . 95 /SISATAN ((A* 77 ./ 80 .- 29 .)si ./15.) 

* +2 . 20/S I SATAN ( (A*7 7 . /80 . -59 • ) s 1 . /40 . ) 

s +1 .55/S I SATAN (- (A*77./80.-7A.)*l ./30.) 

* +2 . 30/S I SATAN (- (A*77 . /80 . - 1 39 • ) * 1 /20 . ) 

s -. 2635 -. 0199) * 2 . 17/2. 1 +. 0199 


CD0N5 = 0.32/SI*ATAN((A*80./85.+8.)*l./30.) 

* + 0 . 60 /SISATAN (- (A* 8 o./ 85 .+ 3 . 5 ) si ./6.5) 

* -0 . 30/S I SATAN (- (A*80 . /85 .+7 • ) * 1 . /90 . ) 

S -0 . 20/S I SATAN (- (A*80 . /85 . -4 .) *1 ./15.) 

S + 1 .95/SISATAN ( (A*8o./85.-28.)*i ./ 15 .) 

s +2 .25/SI SATAN ( (A*80 . /85 . -68 . ) s 1 . /40 . ) 

S +1 . 664/S I*ATAN (- (AS 80 ./ 85 .- 90 .) *1 ./ 30 .) 

S +2 . 35/S I SATAN (- (A*80 . /85 . - HO .) * 1 /20 . ) 

s -.246 


* 


CDON = 0.50/S I SATAN ((A+5.)*l./30.) 

* +0.60/S I *ATAN (- (A+0.) * 1 ,/ 6 .) 

S -0 . 25/S I *ATAN (- (A+3.) *1 ./90.) 

* -0.15/S I SATAN (- (A-4.) *1 ./40.) 

ft + 1 . 85 /S I SATAN ( (A- 30 .) * 1 ./ 28 .) 

* +2 .30/S I SATAN ( (A-60 .) *1 ./40 .) 

* +1.15/SI*ATAN(-(A-85.)/O./30.) 

s +2. 30/S I *ATAN (- (A- 140.) *1/20.) 

* -.2425 


ititititititisititititititititititicitit 


it ic * * j'c j’c j'c j'c * j'c j'c j'c ft j'c jY ic it it ic is is S jY S ft * A j’c * * j'c jY j'c is ic 


CM0N3 = 0.26/SI SATAN (- (A-5. )*l./60.) 
s -0 . 39/S I SATAN ( (A- 1 . ) * 1 . /30 . ) 

* +0.80/SISATAN ((A-5.) *1 ./45.) 

* +0 . 80/S I SATAN (- ( A- 1 0 . ) * 1 . /65 . ) 

* +1 . 80/S I SATAN ( (A-49 .) *1 ./40 .) 

* +2 . 80/S I SATAN (- (A- 69 • ) * 1 /2 3 • ) 

* -0 . 45/S I SATAN (- (A-79 • ) * 1 • /2 • ) 

* -.138 


C 
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CM0NZ3 

= 0.26/SI*ATAN(-(A-5.)*l./60.) 

* 

-0.39/SI*ATAN((A-l . ) * 1 ./14.0) 

sV 

+ 0 . 85 /S 1 *ATAN ( (A~5 •) *1 .7^2.) 

A 

+0 . 80/S 1 *ATAN (- (A-50 . ) * 1 . /60 . ) 

A 

+ 1 .80/S|s':ATAN((A-70.)*i ./5 1 * •) 

>V 

+2 .40/S 1 *ATAN (- (A-69 0*1/25.) 

»’c 

-0 . 45/S 1 *ATAN (- (A-790 *1 ./2.) 

* 

- . 1 38- .02 

CMON53 

=CMON56 

CM0X03 

= CM006 

CMOX23 

= CM0X26 

CMOX53 

= CMOX56 

CM0X3 

= 0 . 26/S 1 *ATAN (- ( A-5 .) * 1 . / 1 0 .) 


-o. 39 /si*atan((a-i 0*1 ./8.) 

* 

+0 .75/SI *ATAN ( (A-5 . ) *1 . / 1 3 0 

>*c 

+ 0 . 70 /S 1 *ATAN (- (A- 10 .) *1 ./ 65 O 

s't 

+ 1 . 20/S l*ATAN((A-49 0*1 -/ 15 0 

5’C 

+2 . 10/S 1 *ATAN (- (A-690 *1/150 

* 

-0 . 45/S 1 *ATAN (- (A-77 0 *1 ./2 0 


-.408+. 01 

CM0X8 

= 0.26/S 1 *ATAN (- (A-50 *1 - / 15 0 

* 

-0 . 33/S 1 *ATAN ( (A - 1 0*1 ./7 0 

s': 

+0. 72/S l*ATAN ((A-5 0*1 ./150 

s't 

+ 0 . 70 /S 1 *AT AN (- (a-350 *1 ./ 750 

s'e 

+ 1 . 13 /SI *ATAN ((A -5 1 0*1 ./I 1 0 

ft 

+2 .08/S 1 *ATAN (- (A-670 *1/170 

3*C 

-0 . 45/S 1 *ATAN (- (A- 78 .) * 1 ./4.) 

* 

-.440 

CMOX 58 

= 0 . 26/S 1 *ATAN (- (A~5 0 *1 . / 10 .) 

* 

"0 . 39/S 1 *ATAN ( (A-l .) *1 . / 1 2 0 


+0. 70/S l*ATAN ((A-5.) *1 . / 15 0 

sV 

+0 . 85 /S 1 *ATAN (- (A-20.) *1 ./40.) 

3*C 

+ 1 . 45/S 1 *ATAN ( (A-52 .) * 1 . / 1 5 0 

* 

+2 . 20/S 1 *ATAN (- (A-69 0*1/150 

sV 

-0.45/SI*ATAN (- (A-77 0*1 ./3- 1 *) 

j'c 

-. 318-. 02 

CMOX28 

= 0 . 36 /SI*ATAN (- (A- 5 .) *1 ./350 

i'c 

-0. 29/S l*ATAN( (A-l 0*1 ./30.) 


+ 1 . 00/S 1 *ATAN ( (A- 1 5 . ) * 1 . /90 .) 

j’c 

+0.75/SI*ATAN (- (A-48.) *1 ./1 10. 

ft 

+0. 90/S 1 *ATAN ((A-52.) *1 ./90 

3*C 

+2 . 1 0/S 1 *ATAN (- (A-69 0*1/1 7 0 

* 

-O. 45 /S 1 *ATAN (- (A-77 0 *1 ./30 

sV 

-.38-. 007 

CM0X08 

= O. 36 /SI *ATAN (- (A-5.) *1 ./350 

* 

-0. 29/S l*ATAN ((A-l 0*1 ./30.) 

* 

+ 1 . 00/S 1 *ATAN ( (A- 1 5 0 * 1 • /90 0 

* 

+0 .80/SI *ATAN (- (A-48 . ) * 1 . /90 . ) 

Vc 

. +0 . 90/S 1 *ATAN ( (A-52 0 * 1 . /9 0 

* 

+2 .10/SI *ATAN (- (A-69 0 * 1 / 1 7 0 

* 

-0 . 45/S 1 *ATAN (- (A-77.) *1 ./30 

s': 

-. 38 -. 007 

CMON58 

= 0. 36 /S l*ATAN(- (A-5 0*1-/350 

* 

-0.29/SI *ATAN ( (A-l .) *1 ./30.) 

sV 

+ 1 . 30/S 1 *ATAN ( (A- 1 5 0 * 1 . /95 0 

3*C 

+0 i 80/S 1 *ATAN (- (A-47 0 * 1 • /35 0 

s': 

+ 1 ,00/S 1 *ATAN ( (A-530 *1 -/ 100 
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* 

* 


+2 . 15/S I *ATAN (- (A- 69 .) * 1 / 18 .) 
-0.45/S I *ATAN (- (A- 78.) *1 ./2 .) 
-.3&--007 


C 


C 


C 

C 


C 


C 


C 


c 


c 


CM 0 NZ 8 = O. 36 /S I *ATAN (- (A-5.) * 1 ./35-) 

* -O. 29 /S I *ATAN ( (A- 1 .) * 1 ./ 30 .) 

* +1 .25/S I *ATAN ( (A- 15.) *1 - /95 - ) 

* +0.80/S I *ATAN (- (A-47 .) *1 ./28.) 

* +1 .00/SI*ATAN ( (A-53 - ) *1 - / 1 0 . ) 

* + 2 . 25/S I *ATAN (- (A- 69 .) * 1 / 18 .) 

* -0.45/S I *ATAN (- (A- 78 .) * 1 ./ 2 .) 

* -.31". 007 

CM0N8 = 0.26/S I *ATAN (- (A-5.) *1 ./40.) 

* -0.45/S I *ATAN ( (A— 4.) *1 ./30 .) 

* +0 . 70/S I *ATAN ( (A-2 .) *1 ./40.) 

* +0 .80/S I *ATAN (- (A— 37 * ) *1 . /25 ♦ ) 

* +1 . 90/S I *ATAN ( (A-52 . ) * 1 . /25 . ) 

* + 2 . 70 /S I *ATAN (- (A- 69 .) * 1 / 20 .) 

* -0.45/SI*ATAN(-(A-79.)*l ./2 .) 

* -. 188 -. 01 


CL 0 X 9 = 0 . 86 /SI AATAN (- (A+ 5 .) * 1 ./ 100 .) 

* +2 ,59/S I *ATAN ( (A~3*) *1 ./7 •) 

* +1 .60/SI*ATAN ( (A-20 .) * 1 . /22 . ) 

* +3.41/S I *ATAN (- (A-57.) *1 ./30.) 

* +.41/S I J-ATAN (- (A-70.) *1 . /20.) 

* -.65 

CL0N9 = 1 .06/SI*ATAN (- (A+5.) *1 ./ 100 .) 

* +1 . 79/S I *ATAN ( (A-5 .) * 1 ./7 .) 

* +2 .50/S I *ATAN ( (A- 1 3 • ) * 1 */22 .) 

* +2.71/SI*ATAN (- (A-59-) *1 ./50.) 

* +1 .21/SI*ATAN (- (A-70.) *1 ./20 .) 

* - . 80 

CM0X9 = 0.26/S I *ATAN (- (A ~5 • ) *1 . / 1 0 . ) 

* -0 . 39/S I *ATAN ( (A-0 . ) * 1 . / 1 0 . ) 

* +1 .00/S I *ATAN ( (A— 3 * ) * 1 ./25-) 

* +0.70/SI*ATAN (‘ (A-7.) *1 ./25.) 

* + 1 . 30 /S I *ATAN ( (A-50 . ) * 1 . / 1 6 . ) 

* + 2 . 10 /S I *ATAN (- (A- 69 .) * 1 / 15 .) 

* -0.45/S I *ATAN (- (A- 76 .) * 1 ./3 -5) 

* -.433 

CMOX59 = 0.26/S I *ATAN (- (A-5.)*l ./4 . ) 

* -0 .39/SI *ATAN ( (A+2 . ) * 1 . /5 • ) 

* +1 . 00/S I *ATAN ( (A- 1 . ) * l . / 15 . ) 

* + 0 . 60 /SI*ATAN (- (A- 18.) * 1 / 13 .) 

* + 1 . 30/S I *ATAN ( (A-5 1 . ) * 1 • / 1 4 . ) 

* + 2 . 15/SI *ATAN (- (A- 69 .) *1/15.) 

* -0 . 45/S I *ATAN (- (A-76.) *1 ./3 .5) 

* -.490 

CMOX29 = 0.26/SI*ATAN (- (A-5.) *1 ./10.) 

* -0 . 39/S I *ATAN ( (A+4 .) * 1 ./7.) 

* +1 .33/SI*ATAN ( (A— 1 .) *1 ./40 .) 

* +0.60/SI*ATAN (- (A- 19 .) *1/17.) 

* +0.90/S I *ATAN ( (A-5 1 •) *1 ./9-) 

* + 2 . 05 /S I *ATAN (- (A- 69 .) *1/15.) 

* -0.45/S I *ATAN (- (A- 76 .) * 1 ./3-5) 

* -.450 

CMOX09 = 0.26/S|5'<ATAN (~ (A~5.) *1 . / 1 0 . ) 

* -0.39/SI *ATAN ( (A+2.) *1 ./10.) 
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c 


c 


c 

c 


c 


c 


c 


c 


c 


* +1.30/SI*ATAN((A-1.)*1. 1/1*0.) 

* +0.60/SI*ATAN (- (A-19-) *1/15 -) 

* + 0 . 90 /S I *ATAN ( (A~51 .) *1 ./9 •) 

* +2 . 1 1 /S I *ATAN (- (A-69 • ) * 1 / 1 5 • ) 

* -0.45/S I *ATAN (- (A-7&.)*1 . /3 -5) 

* -.490 


CMON59 = 0.26/SI*ATAN(-(A-5.)*l ./I 1 .) 

* -0 . 39/S I * ATAN ( (A+2 . ) * 1 . / 1 0 . ) 

A +1 .50/SI AATAN ((A-l .) *1 . 1/48.) 

A +0 . 60/S I AATAN (- (A- 1 9 . ) * 1 / 1 5 . ) 

* +0. 80/SI AATAN ((A-51.)*l./10.) 

* + 2 . 32 /S I AATAN (- (A- 70 .) * 1 / 20 .) 

A -0.45/S I AATAN (- (A- 78 .) *1 . /3 . 5) 

* -.490 


CM0NZ9 = 0.26/S I AATAN (-(A-5.)Al./7.) 

A -0. 39/S I AATAN ((A+0.)*1 ./10.) 

A +1 .60/SI AATAN ((A-5.)*l ./50.) 

* +0 . 60/S I *ATAN (- (A— 32 . ) *1/1 1 .) 

A +0.80/Sl*ATAN((A-51 .) *1 ./I 1) 

* +2. 23/S I *ATAN (- (A- 69 .) * 1 / 19 .) 

* -0.45/S I *ATAN (- (A- 78 .) * 1 . /3 . 5) 

A -.410 


CM0N9 = 0.l6/SI*ATAN(-(A-5.)*l./40.) 

* -0 . 39/S I AATAN ( ( A- 3 . ) * 1 . / 8 . ) 

A +1 . 20/S I AATAN ((A- 15.) * 1 ./ 1 20 .) 

* +0 . 70/S I *AT AN (- (A— 25 • ) *1 ./50.) 

A + 2 . 00 /S I *ATAN ( (A- 52 .) A1 ./75.) 

* +2. 00/S I AATAN (- (A-69.) *1/20.) 

A -0.42/Sl*ATAN (- (A- 79 .) *1 ./It.) 

A -.068 

CLQ = 0.26/SI *ATAN (- (A~5.) *1 ./10.) 

* -2 . 39/S I AATAN ( (A- 6 . ) * 1 . /3 . ) 

A +2 . 40/S I AATAN ( (A- 1 5 . 5) * 1 . /3 . ) 

A +2 . 00/S I *ATAN (- (A- 20 . ) A 1 . /5 . ) 

* +A. 30/S I AATAN ( (A-37.) *1 ./4. 5) 

A +2 . 20/S I *ATAN (- (A-45-) *1/15-) 

A +2. 20/S I *ATAN (- (A-80.) *1/15-) 

A - 0 . 45/S I *ATAN (- (A- 76 .) * 1 ,/3 .5) 

A +4.2 

CLAD = 1 .32/P I *ATAN (50 . *P I / 1 80 . * (-A+5 . ) ) + 1 .8 
A -0.75*ATAN(l.*(A-A5.)/2) 

CMQ =-0 . 82/P I AATAN (20 . *P I / 1 80 . * (-A+5 . ) ) ~5 • 8 

* +2 . 00*AT AN ( (- A+32) / 6 . ) 

it +4.55AATAN ((A-43.) *3-5) 

A - 3 . 5 OAATAN ( (A-57 • ) /5.) 

CMAD =-0 .02/P I AATAN (50 . *P I / 1 80. * (-A+1 .)) -0.9 
A +0.50AATAN ((A- 6 .) *5) 

A -0.80*ATAN ( (A— 18.) /2) 

A +0 . 90AATAN ( (A — 45 • ) /2) 

A 1 =CL0X6 
A2 =CL0N6 

A 1 1= (A1-A2) / (10.5- (-24) ) * (DH- (-24) )+A2 
A3 =CL0X9 
A4 =CL0N9 

A 1 2 = (A3"A4) / ( 1 0 . 5 - (-24) ) A (DH- (-24) ) +A4 
CLO = All + (A 1 2-A 1 1) * (MACH- . 6 ) /0.3 
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CLO=A 1 1 IF MACH NUMBER =.6 

CLO=A 1 2 IF MACH NUMBER = .9 


IF (DH.GE.5-) THEN 

BIO =(CMOX3"CMOX53) *(DH-5.)/5.5+CMOX53 
B 1 1 =(CMOX6-CMOX56)*(DH-5.)/5.5+CMOX56 
B 1 2 = (CMOX8-CMOX58) * (DH-5 . ) /5 . 5+CMOX58 
B 1 3 = (CMOX9-CMOX59) * (DH-5 . ) /5 ■ 5+CMOX59 
GO TO 200 

ELSE 

GO TO 5 
END IF 
5 CONTINUE 

IF (DH.GE.2.) THEN 

BIO =(CMOX53-CMOX23) * (DH-2 . ) /3 .+CMOX23 
B 1 1 =(CMOX56-CMOX26)*(DH-2.)/3.+CMOX26 
B 1 2 = (CMOX58-CMOX28) * (DH-2 .) /3 .+CMOX28 
B 1 1 =(CMOX59-CMOX29)*(DH-2.)/3.+CMOX29 
GO TO 200 

ELSE 

GO TO 8 
END IF 
8 CONTINUE 

IF (DH.GE.O.) THEN 

BIO = (CM0X23-CM0X03) *DH/ 2 .+CM0X03 
B1 1 = (CM0X26-CM006) *DH/2 .+CM006 
B 1 2 = (CM0X28-CM0X08) *DH/2 .+CM0X08 
B 1 3 = (CM0X29-CM0X09) *DH/2 .+CM0X09 
GO TO 200 

ELSE 

GO TO 10 
END IF 
10 CONTINUE 

IF (DH.GE.-5.) THEN 

B 1 0 = (CM0X03-CM0N53) * (DH+5 . ) /5 .+CM0N53 
B 1 1 = (CMOO 6 -CMON 56 ) * (DH+ 5 .) / 5 .+CMON 56 
B 1 2 = (CM0X08-CM0N58) * (DH+5 . ) /5 -+CM0N58 
B 1 3 = (CM0X09-CM0N59) * (DH+5 • ) /5 -+CM0N59 
GO TO 200 

ELSE 

GO TO 20 
END IF 
20 CONTINUE 

IF (DH.GE.-12.5) THEN 

B 1 0 = (CMON53-CMONZ3) * (DH+ 1 2 . 5) /7 -5+CMONZ3 
B 1 1 = (CMON56-CMOZ6) * (DH+1 2 .5) /7-5+CMOZ6 
B12 =(CMON58-CMONZ8) * (DH+12.5)/7-5+CMONZ8 
B 1 3 = (CMON59-CMONZ9) * (DH+1 2 .5) /7 - 5+CMONZ9 
GO TO 200 

ELSE 

GO TO 30 
END IF 
30 CONTINUE 

IF (DH.GE.-2lt.) THEN 

BIO = (CMONZ3-CMON3) * (DH+21+ .) / 1 1 .5+CMON3 
B 1 1 = (CMOZ 6 -CMON 6 ) * (DH+21* .) /1 1 . 5 +CMON 6 
B 12 =(CMONZ 8 -CMON 8 )*(DH+ 2 lt.)/l 1 . 5 +CMON 8 
B 1 3 = (CMONZ9-CMON9) * (DH+21* .) / 1 1 .5+CMON9 

ELSE 

GO TO 200 
END IF 

200 CONTINUE 

C 

IF (MACH. LE. 0.6) THEN 

CMO = B 10 + (B 1 1 -BIO) * (MACH- . 3 ) /O . 3 
GO TO 250 
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210 


ELSE 

GO TO 210 
END IF 
CONTINUE 

IF (MACH. LE. 0.8) THEN 
CMO = Bll + (B 1 2-B 1 1) * (MACH- .6) /0 . 2 
GO TO 250 
ELSE 

GO TO 220 
END IF 
220 CONTINUE 

IF (MACH.LT.l.) THEN 

CMO = B12 +(B13-B12)*(MACH-.8)/0.1 

ELSE 

GO TO 250 
END IF 

250 CONTINUE 


CMO = BIO IF MACH NUMBER =.3 
CMO = Bll IF MACH NUMBER = .6 
CMO = B12 IF MACH NUMBER =.9 


IF (DH.GE.O.) THEN 

CDO =(CD0X-CD0Z)*(DH+0.)/10.5+CD0Z 
C PRINT *,'DH IS LAGER THAN O. 1 
GO TO 300 

ELSE 

GO TO 50 
END IF 
50 CONTINUE 

IF (DH.GE.-5.) THEN 

CDO = (CD0Z-CD0N5) * (DH+5-) /5-+CD0N5 
C PRINT *,'DH IS LAGER THAN ~5- ' 

GO TO 300 

ELSE 

GO TO 60 
END IF 
60 CONTINUE 

IF (DH.GE.-2A.) THEN 

CDO = (CD0N5-CD0N) * (DH+2A.) /19.+CD0N 
END IF 

300 CONTINUE 

c Mmmmmmmmmmmmmmmmmmmmm 

C end of interpolations 

RETURN 

END 
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C 


c— 



— 




-“== 

c 

PURPOSE: 

THIS 

SUBROUTINE WILL BE CALLED 

TO CALCULATE 12 AERODY- 

== 

c 


NAMIC 

COEFFICIENTS WHICH ARE THE OUTPUTS 

OF THIS SUB- 

== 

c 


ROUTINE (LATERAL COEFFICIENTS) 



= = 

c 

INPUTS: 






SC s 

c 


ALPDEG 

ANGLE OF ATTACK 

(DEG) 


as 

c 


BETDEG 

SIDESLIP ANGLE 

(DEG) 


== 

c 


MACH 


MACH NUMBER 



— 

c 


DA 


AILERON DEFLECTION 

(DEG) 


== 

c 


DR 


RUDDER DEFLECTION 

(DEG) . 


== 

c 


DH 


STABILATOR DEFLECTION 

(DEG) 


== 

c 


P 


AIRCRAFT X BODY AXIS 

ROLL RATE 

(RAD/SEC) 

== 

c 


R 


AIRCRAFT Z_BODY AXIS 

YAW RATE 

(RAD/SEC) 

== 

c 

OUTPUTS: 






SS 

c 


CRO 

, CYE 

1 , CYP , CYR 



== 

c 


CIO 

, Cl E 

1 ,C1P ,C1R 



sees 

c 


CNO 

, CNE 

1 , CNP , CNR 



SB 


C == 

C THOSE COEFFICIENTS ARE USED IN AEROLAT IN THE FOLLOWING FORMULAS == 
C 

C CY = CYO+CYB*BETDEG+B02VT*(CYP*P+CYR*R) 

C Cl = C10+C1B*BETDEG+B02VT*(C1P*P+C1R*R) 

CN = CN0+CNB*BETDEG+B02VT* (CNP*P+CNR*R) 

B02VT = BREF/ (2*V) 

CY SIDE FORCE COEFFICIENT ALONG Y_WIND AXIS 

Cl ROLL MOMENT COEFFICIENT ABOUT X_WIND AXIS 

CN YAW MOMENT COEFFICIENT ABOUT Z WIND AXIS 

BREF REFERENCE WINGSPAN 

V AIRCRAFT TOTAL AIRSPEED 

ss 

AUTHER: J I CHANG CAO 

GRADUATE RESEARCH ASSISTANT 
GEORGIA INSTITUTE OF TECHNOLOGY 
ATLANTA, GA30332 


SUBROUTINE NAER0C2 (CYO , C 1 0 , CNO ,CYB ,C1B ,CNB , 

* CYR , C 1 R , CNR , CYP ,C1P ,CNP , 

=F ( 

* ALPDEG,BETDEG,MACH,DA,DR,DH,P,R, ALT) 
IMPLICIT REAL (C) 

REAL ALFA,BETA,MACH,DA,OR,DH,P,R, ALT, ALPDEG.BETDEG, A, 

* DAN , DAX, DRN, DRX, DHN , DHX, P I , HC 1 P , LC 1 P 


CONVERSION OF ALFA AND BETA 


IF 

(DA.LT.-25.) 

DA=-25. 

IF 

(DA.GT.25.) 

DA=25 • 

IF 

(DR.LT.-30.) 

DR=-30. 

IF 

(DR.GT. 30 .) 

DR= 30. 

IF 

(DH.LT.-2L.) 

DH=-2L. 

IF 

(DH.GT.10.5) 

DH=10.5 

A 

= ALPDEG 


PI 

=ACOS (-1 

.) 

SI 

=2.75 



EXTREMAL VALUES OF CONTROL SURFACES DEFLECTION 
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DAN = -25. 
DAX = -DAN 
DRN = -30. 
DRX = -DRN 
DHN = -21*. 
DHX = 10.5 


COMPUTATION OF COEFFICIENTS 


CYOXXBO = .01216/S I AATAN (A*3* A •) 

+ - 03247/S I AATAN (- (A- 1 3 .) /I* .) 

+ .OO 89 I/SI *ATAN ( (A-29.5) *2-/3-) 
+ .03058/S I *ATAN ( (1*6 . - A) *2 . /5 • ) 

+ .02759/S I *AT AN ( (75 • -A) *3 • AO .) 
+.03477 

CYOXXB2 = .06/PhVATAN (A/3.) 

+ . 09/PI *ATAN ( ( A — 3 1 .) *5./8.) 

+ .06/P I AATAN (- (A-46.) *3./ 10.) 

+ .03/PI*ATAN ((A- 63 .) *3./7.) 

+ .09/PI *ATAN ( (75 - ~ A) *3./5.) 

+ .04/PI AATAN ( (A- 85 .) *4./5.) 
-0.285 

CYOXNBO = 0.029624*ATAN ((A-20.) *4./25.) 
-0 . 0020868 AATAN ( ( 1 2 . -A) *5 . ) 
+6.99*EXP (- (A+17 . 6 ) ) 

-0.075535 


CY0XNB2 = .232/PIAATAN (1 10*PI/180.*(A-16.)) -.394 
CYONXBO * -CYOXNBO 


CY0NXB2 = .047/PIAATAN (A/2.) 

+.021 /PI AATAN ( (17- —A) *5 . /4 . ) 
+ .037/P I *ATAN ( (A-32 . ) *5 • /3 • ) 
+ . 06 /P I *ATAN ( ( 76 . -A) *5-/4.) 
+.043/PIAATAN (A- 85 .) 

-.248 


* 

A 


CYONNBO = -CYOXXBO 

CY0NNB2 = .219544/P I AATAN ( (A-21 .) *5./56.) 
+ . 0636375/P I *ATAN (- (A-74.) /4.) 

+ . 0709 125/P I AATAN ((A- 85 .) A 3 . / 10 .) 
- . 36444 

C10XXB0 = .04226/SI AATAN (- (A-20.) A2./7.) 

+ .OO 83 I/SIAATAN (- (A- 53 .) Al*./ 7 .) 

+ . 00997/S I *ATAN ( (A-65 .) Ai* ./ 5 .) 

+ .010 1/S I AATAN (- (A- 77 . 5 ) A 8 ./I 5 .) 
-.002/PI AATAN (- (A— 8 . ) A 10.) 

+.0286 


C10XXB2 = .0277 1/S I AATAN (— (A — 2 . ) /5 • ) 

> +. 06763/S I AATAN (- (A- 19.) /4.) 

+ .006/S I AATAN ( (A— 22 . ) *2.) 

+.00081 


* 

A 


C10XNB0 = .0041/S I *ATAN ( (A-4.) *2./3.) 

+ .003/SIAATAN (- (A-20.)/3.) 

+ .01 1/SI AATAN (- (A-59.)*2./5.) 
-.00144/S I AATAN (- (A+l .)a 8 .) 

+.01793 

C10XNB2 = .085/SI*ATAN(-(A-15.)a8./92.)-.0065 
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c 

c 


c 

c 


c 


c 

c 


c 


c 

c 


c 

c 


C10NXB0 = -C10XNB0 


C10NXB2 = .0395/SI*ATAN (- (A-5-) *3./13.) 

* + .0295/S I *ATAN ( (A-25 .) *3 • H ■ ) 

* +. 01 26 /S I *ATAN (- (A- 38 .) *4 ./3 •) 

* +.01 l4/SI*ATAN((A-42.)*2.) 

* +. 0082/S I *ATAN {- (A-51 •) /2 .) 

* +.0 1 32/S I *ATAN ( (A— 70 .) *2./5.) 

* -.0479 

C10NNB0 = -C10XXB0 

C 10 NNB 2 = .0344/SI*ATAN(-(A-5.) *3-/13.) 

* +. 037/S I *ATAN ( (A-24 .) *2 ./5 •) 

* +. 01 1 /S I *ATAN (- (A- 38 .) * 2 .) 

* +.01 2/S I *ATAN ( (A-42 .) *7-/4.) 

* +.01 1/S I *ATAN (- (A— 52 .) *3-/8.) 

* +.0176/SI*ATAN ( (A— 73 •) *3 ./l 3 .) 

* “*0553 

C 10 XX 0 = (.04226/SI*ATAN(-(A-8.5)*2./7.) 

* +.01031/S I *ATAN (- (A- 35 - ) *4-/7.) 

* +.00997/S I *ATAN ( (A - 65 .) *4 ./ 5 *) 

*. + .0 1 31/S I *ATAN (- (A- 77 . 5 ) * 8 . / 1 5 . ) 

* - .002/P I *ATAN (- (A-8.)*10.) 

* +.0286) /l. 8 

CL0XX2 = .085/S l*ATAN (- (A-15 •) *8./92 .) - .01 98 

C10XN0 =( .004 1/S I *ATAN ( (A— 8 .) *2-/3-) 

* +.0 1 2/S I *ATAN (- (A- 1 1 . ) / 3 • ) 

* +.01 0/S I *ATAN (- (A-15 •) /3.) 

* +.0 1 2/S I *ATAN ( - (A~35 •) /3 • ) 

* +. 005/S I *ATAN (- (A- 100.)*: ./5.) 

* -.00144/SI *ATAN (- (A-2.)*8.) 

* +.01793) /2.~. 0032 

C10XN2 = .048/S I *ATAN (- (A-l 8 .) *8./92 .) 

* +. 030/S I *ATAN (- (A-23-)*8./92.) 

* +. 045/S I *ATAN ( (A-80.) *8./92 .) 

* -.0105 

C10NX0 = -C10XN0 

C 10 NX 2 = 0.0295/SI*ATAN(-(A-25.)- ; 3-/13.) 

* + 0 .0295/S I *ATAN ( (A-42 . 5) *3 • /7 • ) 

* +. 0086/S I *ATAN (- (A-15.) *4. /3.) 

* +.001 4/S I *ATAN ( (A-42.) *2.) 

* +.01 62/S I *ATAN (- (A-51 -)/2.) 

* + . 00 1 2 /S I *ATAN ( (A-70 .) *2 . /5 •) 

* -.0479+- 00 13 

C10NN0 = -C10XX0 


C10NN2 =- (0.0544/S I *ATAN (- (A~7 . ) *3 ./ 1 3 .) 

* +0 . 087/S I *ATAN ( (A— 17.) *2./7.) 

* +.008/SI*ATAN (- (A-25.) *2.) 

* +. 0 1 9/S I *ATAN ( ( A-28 . ) *7 . /4 . ) 

* +. 051/S I *ATAN (- (A-42.) *3-/8.) 

* +. 0096/S I *ATAN ((A-55.) *3-/13.) 

* -.0553) /2. 6-0. 074 

CNXXXBO = .00952/S I *ATAN ( (A— 1 4 . ) *3 • /8 .) 

* +.0 1056/S I *AT AN ( (A-47 •) *2 ./3 .) 

* +. 01 395/S I *ATAN((A-67.) *9-/14.) 
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ft 

ft 

C 

CNXXXB2 

ft 

ft 

ft 

ft 

C 

CNXXNBO 

ft 

ft 

C 

CNXXNB2 

ft 

ft 

ft 

ft 

C 

CNXNXBO 

ft 


C 

CNXNXB2 

ft 

* 

ft 

ft 

C 

CNXNNBO 

* 

ft 

ft 

C 

CNXNNB2 


C 

C 


CNNXXBO 

CNNXXB2 

ft 

ft 


C 

C 

C 

C 



+ . 00899/S I *AT AN (- (A -81 .) *3-/8.) 
-.01862 

= . 077832 /s I *ATAN (- (A- 27 . 5 ) *4./45.) 
+ . 0744/S I *ATAN ((A-58.5) *4./25*) 

+ . 02885/S I *ATAN (- (A-79 •) *7 - /20 .) 

+ . 006/S I *ATAN (- (A-43 •) ) 

-.022 

= .00953/SI*ATAN((A-15-)*3./8.) 

+ .004 1 1 /S I *ATAN (- (A-46 . 5) *8 . / 1 5 • ) 
+ . 02222/S I *ATAN((A-71 .5) *4-/25-) 
-.01781 

= . 080916 /S I *ATAN (- (A- 30 .) *3-/42.) 
+ . 056/S I *ATAN ( (A-62 . ) /4 . ) 

+ . 02085/S I *ATAN (- (A-79 -) /5 -) 
-.005/SI*ATAN (- (A-82.) *2.) 

-.02221 

= .02026/S I *ATAN (- (A- 19.) *3-/ 1*4.) 

+ . 022/S I *ATAN (- (A-i+9 .) *2./7.) 

+ • 03393/S I *ATAN ( (A-73.) /7 • ) 

+ . 002/S I *ATAN (- (A-14.) *2.) 

+ . 003 /S I *ATAN ( (A- 76 .) *2 .) 

+.02769 

= . 06482/S I *ATAN (- (A-20 . ) *3 . / 1 4 . ) 

+ .Ol*937/SI*ATAN (- (A-41 .) *5- /I 1 *-) 

+ -053/S I *ATAN ( (A-60.) / 6 .) 

+ . 005/S I *ATAN (- (A-8.) *4.) 

+.0211 

= .02037/S I *ATAN (- (A- 1 7 • 5) * 8 . A3 • ) 
+ . 00389 /s 1 *ATAN (- (A-54 . 5) * 6 . / 1 7 . ) 
+ . 0 1 623 /S I *ATAN ( (A- 69 . 5 ) *5-/23-) 

+ . 002/S I *ATAN (- (A- 13.) *2.) 

+.02711 

= .1 1857/S I *ATAN (- (A-26.5) *4./45.) 

+ . 07065/S I *ATAN ((A-60. 5) A./25.) 

+ .025 l 8 /S I *ATAN (- (A-80.) *3-/10.) 

+ .005/SI*ATAN (- (A- 18 .) * 2 .) 

+.020265 

= -CNXNXBO 

= .06 1 32/S I *ATAN (- (A— 3 1 - ) / 1 0 . ) 

+ .0552 1 /S I *ATAN ( (A— 55 • ) *2./9-) 

+ . 04659/S I *ATAN (- (A- 77 • 5) A • / 1 5 • ) 
-.03235 

= -CNXNNBO 

= .O 63768 /S I *ATAN (- (A — 32 .) / 1 0 .) 

+ .04788/S I *ATAN ( (A-57 • 5) * 6 . /25 • ) 

+ . 03829/S I *ATAN (- (A-77-5) *4./ 15.) 
-.03288 

= -CNXXXBO 

= • 1 1977/S I *ATAN (- (A-28.5) *3-5/51.) 
+ . 04303/S I *ATAN ((A-58.5) *4. /1 3.) 

+ . 02532/S I *ATAN (" (A-74.) *2./5.) 
+.01184 


C 
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c 


CNNNNBO = -CNXXNBO 


C 


C 


C 


C 


C 


C 


C 


C 


C 


CNNNNB2 = .120543/SlftATAN(-(A-28.)*3./40.) 

* +.05707/SI*ATAN((A-55.5)*W25.) 

* +. 03650 /S I *ATAN(- (A- 78 . 5 ) *4 . / 1 5 * ) 

* +.01202 

CNXXXO =( .011*52/SI*ATAN((A-ll.)*3./8.) 

* -.005/S I *ATAN((A-22.) *3-/8.) 

* -t-.Ol 156/S I *ATAN ( (A-i*l .)*2./3.) 

* +.01 205/S I *AT AN ( (A-67 • ) *9 • / 1 4 . ) 

* +. 00769/S I *ATAN (- (A- 8 l .) * 3 . / 8 . ) 

ft - .01862) /I .2 1 

CNXXX2 = .067832/S I *ATAN(-(A-jO.) *4. A5) 

* +.O 65 /S I *ATAN (- (A-l4.)*4./45.) 

* - .04/S I *ATAN (- (A-4.)*4./45.) 

* +. 0844/S I ftATAN ( (A-46 .) *4 ./25 •) 

* +. 04085/S I ftATAN <- (A- 100 .) *7 ./20 .) 

* +.006/SI*ATAN(-(A-35.)) 

ft -.0352 

CNXXNO = . 01253/SI *ATAN ((A-12.) *3. / 8 .) 

ft -. 00200/S I *ATAN ( (A-22 .) *3 . / 8 .) 

* +.0041 1/SI*ATAN (- (A-46. 5) *8./15.) 

* +. 02222/S I *ATAN ( (A — 8 1 . 5 ) *4./25.) 

* -.01161 


CNXXN2 = . 1 1091/S I *ATAN (- (A-25-) *3-/42 .) 

* - . 025/S I *ATAN (- (A-8.) *3-/42.) 

* +. 05600/S I *ATAN((A-48.)/4.) 

* +. 03385/S I *ATAN (- (A— 1 00 . ) /5.) 

* - . 005/S I *ATAN (- (A— 82 . ) *2 . ) 

* -.03125 

CNXNXO = .01026/SI *ATAN (- (A- 13 •) *3 • / 1 4 .) 

* - . 009/S I *AT AN (- (A-40.)*2./7.) 

* +.010/SI*ATAN (" (A- 18.) *2./7.) 

* - .002/S I *ATAN (- (A-30.) *2./7.) 

* +. 022/S I *ATAN (- (A-49 .) *2 ./7 .) 

* +. 02543/S I *ATAN ( (A-83 -) /7 •) 

* +. 002/S I *ATAN (- (A-7.) *2.) 

* +. 003/S I ftATAN ( (A — 76 . ) *2 .) 

ft +.02769 

CNXNX2 = .05428/S I ftATAN (- (A- 16.) *3-/14.) 

* +. 05037/S I ftATAN (- (A-36.) *5./ 14.) 

ft -.003/S I *ATAN ((A-45J/6.) 

* +. 060/S I *ATAN ((A-50J/6.) 

* +.005/S I ftATAN (- (A-8.) *4.) 

ft +.0211 


CNXNNO = .01637/SI ftATAN (— (A— 1 0 . ) *8 . /4 3 - ) 

* + . 00559/S I ftATAN (- (A- 1 80 . ) *6 . / 1 7 • ) 

* +. 0 1 62 3/S I ftATAN ( ( A- 1 00 . ) *5 . /2 3 . ) 

ft +. 002/S I ftATAN (-(A- 13.) *2.) 

* +.0271 

CNXNN2 = .05428/S I ftATAN (- (A- 15 •) *3-/14.) 
ft +• 05037/S I ftATAN (- (A-33 • ) *5 . / 1 4 . ) 

* -.003/S I ftATAN ( (A-45-) /6 .) 

ft +. 062/S I *ATAN ( (A-46. )/6.) 

* +. 005/S I ftATAN (- (A-8 .) *4 .) 

ft +.0241 


CNNXXO = -CNXNXO 
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c 

CNNXX 2 = .O 67832 /S I *ATAN (- (A- 32 .) *4 ./45) 

* +. 063/S I *ATAN (- (A-14.) *4./4 5-) 

* -.04/SI*ATAN (- (A-4.) *4./4 5.) 

* +. 0844/S I *ATAN ( (A— 46 . ) *4 ./25-) 

* +. 04085/S I *ATAN (- (A- 90 .) * 7 -/ 20 .) 

* +.006/SI*ATAN(-(A-35.)) 

* -.0352 

c 

CNNXNO = -CNXNNO 
C 

CNNXX 2 = .O 67832 /SI *ATAN (- (A— 35 - )-*4 - 

* +.070/SI*ATAN (- (A-14-.) *4-/45-) 

* - .04/S I *ATAN (- ( A -4 . ) *4./45.) 

* + . 0844/S I *ATAN ( (A-46 . ) *4 . /25 • ) 

* +. 04085/S I *ATAN (- (A-IOO.) *7-/20.) 

* +. 006/S ! *ATAN (- (A— 35 * ) ) 

* -.0392 
C 

CNNNXO = -CNXXXO 
C 

CNNNX2 = .05428/S I *ATAN (- (A — 2 2 - ) *3./ 14.) 

* +. 02037 /S I *ATAN (- (A- 36 .) * 5 . / 1 4 . ) 

* - .003/S I *ATAN ( (A— 45 - ) / 6 . ) 

* +. 042/S I *ATAN ( (A— 47 . ) / 6 .) 

* +. 020/S I *ATAN (- (A- 11 .) * 4 .) 

* +.0181 
C 

CNNNNO = -CNXXNO 
C 

CNNNN2 = CNXNN2 
C 

CYB = .206*E-3 
C 

CYP = . 086/P I *AT AN ( 1 00 . *P I / 1 80 . A) 

* + . 096/P I *ATAN ( 1 00 . *P I / 1 80 . * (- A+2 3 - ) ) 

* +. 22/P I *ATAN (100 . *P I / 1 80 . * (-A+45 •) ) 

* + . 256/P I *ATAN ( 1 00 . *P I / 1 80 . * (A-54 . ) ) - . 047 
C 

CYR = .17/PI *ATAN (100.*P I /1 80 .* (A-4 .) ) 

* +.55/PI*ATAN (100.*PI/l80.*(-A+20.)) 

* + . 54 /P I *AT AN ( 1 00 . *P I / 1 80 . * (A-45 • ) ) 

* + . 26/P I *ATAN ( 1 00 . *P I / 1 80 . * (-A+6 1 . ) ) +.07 
C 

C1B = 1 . E-4* (6.32/P I *ATAN (1000 .*P I /180.* (-A+1 3-) ) +3-26) 
C 

LC1P = . 15/P I *ATAN (100.*P 1/180.* (A- 12) ) - 

* +. 25/P I *AT AN (1000. *P I / 1 80 . * (- A+ 28 . ) ) 

* +.55/PI*ATAN (1000.*P 1/180.* (A-4 1 .)) 

* +• 33/P I *ATAN (100. *P l/l80. * (-A+50.) ) -.341 
C 

HC1P = .28/PI*ATAN(100.*PI/l80.*(A-10.)) 

* + . 25/P I *AT AN ( 1 000 . *P I / 1 80 . * (- A+4 1 . ) ) 

* +• 55/P I *ATAN (1000 . *P I / 1 80 . * (A-4 1 .) ) 

* + . 33/P I *ATAN ( 1 00 . *P I / 1 80 . * (- A+50 . ) ) - . 47 1 

c 

C1P = LC1P + (HClP-LClP)*(MACH-.6)/0.3 
C — — 

C1R = .304/PI*ATAN(100.*PI/l80.*(A-3.)) 

* + . 2 2/P I *ATAN (20 . *P I / 1 80 . * ( 50 . - A) ) 

* - .026*EXP ( (A- 85 .) / 100 .) +.018 
C 

CNB = 1 . E-6* (12.7 /P I *ATAN (1000.*PI/l80.* (-A+13-) ) “1 1 -7) 
C 

CNP = -075/PI*ATAN (50.*PI/l80.* (A-17-) ) 

* + . 04/P I *ATAN ( 1 00 . *P I / 1 80 . * (A-50 . ) ) 
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o o o o o o 


* + . 2/P I *ATAN ( 1 000 . *P I / 1 80 . * (- A+57 • ) ) 

* + . 1 3/P I *ATAN ( 1 000 . ftp I / 1 80 . * (A-62 . ) ) 

* +.09/PI*ATAN(1000.*p|/l80.*(-A+73.)) 

* +. 1/PI*ATAN(1000.*PI/180.*(a-77.)) -.028 
C 

CNR = . 1 6/P I *ATAN ( 100 . *P I /I 80 . * (-A+22 .) ) 

* +. 34/P I *AT AN ( 100 . *P I / 1 80 . * (A-57 •) ) 

* -.1*EXP((A-78.)/10.)-.09 

c 

c mMmmmmmmmmmmmmmmmmmmmmm# 

A 1 =BETDEG/20ftCY0XXB2+ (20-BETDEG) /20*CY0XXB0 
A2=BETDEG/20*CY0NXB2+ (20-BETDEG) /20ftCY0NXB0 
A 1 1 = (A 1 -A2) / (25- (-25) ) * (DA- (-25) ) +A2 
A3=BETDEG/20ftCY0XNB2+ (20-BETDEG) /20ftCY0XNB0 
A4=BETDEG/20*CY0NNB2+ (20-BETDEG) /20*CY0NNB0 
A 1 2= (A3-A4) / (25- (-25) ) * (DA- (-25) )+A4 
A 1 3= (A 1 1 -A 1 2) / (30- (-30) ) * (DR- (-30) ) +A 1 2 


CY0=A 1 3 


D1=BETDEG/20ftCNXXXB2+ (20-BETDEG) /20ftCNXXXB0 
D2=BETDEG/20*CNXXNB2+ (20-BETDEG) /20ftCNXXNB0 
D3=BETDEG/20*CNXNXB2+ (20-BETDEG) /20*CNXNXBO 
D4=BETDEG/20*CNXNNB2+ (20-BETDEG) /20ftCNXNNB0 
D5=BETDEG/20*CNNXXB2+ (20-BETDEG) /20ftCNNXXB0 
D6=BETDEG/20ftCNNXNB2+ (20-BETDEG) /20*CNNXNB0 
D7=BETDEG/20ftCNNNXB2+ (20-BETDEG) /20ftCNNNXB0 
D8=BETDEG/20ftCNNNNB2+ (20-BETDEG) /20ftCNNNNB0 
D 1 1= (D 1 -D2) / (10 .5“ (-24) ) * (DH- (-24))+D2 
D 1 3== (D3-DA) / (10.5- (-24) ) * (DH- (-24) )+D4 
D15= (D5-D6) / (10.5- (-24) ) * (DH- (-24) )+D6 
D 1 7= (D7-D8) / (10-5- (-24) ) * (DH- (-24) ) +D8 
D2 1 » (D 1 1 -D 1 3) / (30- (-30) ) * (DR- (- 30) ) +D 1 3 
D22*= (D 1 5 -D 1 7) / (30- (-30) ) * (DR- (-30) ) +D 1 7 
D3 1- (D2 1 -D22) / (25- (-25) ) * (DA- (-25) ) +D22 


CNO=D 3 1 IF MACH NUMBER =.6 


DD1=BETDEG/20*CNXXX2+ (20-BETDEG) /20ftCNXXX0 
DD2=BETDEG/20ftCNXXN2+ (20-BETDEG) /20ftCNXXN0 
DD3=BETDEG/20*CNXNX2+ (20-BETDEG) /20*CNXNX0 
DD4=BETDEG/20*CNXNN2+ (20-BETDEG) /20ftCNXNN0 
DD5=BETDEG/20*CNNXX2+ (20-BETDEG) /20ftCNNXX0 
DD6=BETDEG/20ftCNNXN2+ (20-BETDEG) /20*CNNXN0 
DD7=BETDEG/20*CNNNX2+ (20-BETDEG) /20ftCNNNX0 
DD8=BETDEG/20*CNNNN2+ (20-BETDEG) /20ftCNNNN0 
DD 1 1= (DD1-DD2) / (10.5" (-24) ) ft (DH- (-24) )+DD2 
DD 1 3 = (DD3-DD4) / (10 . 5 - ("24) ) * (DH- (-24) ) +DD4 
DD 15= (DD5-DD6) / (10 .5- (-24) ) ft (DH- (-24) ) +DD6 
DD 1 7= (DD7-DD8) / (10.5- (-24) ) ft (DH- (-24) )+DD8 
DD21= (DD1 1-DD13) / (30- (-30) ) * (DR- (~30) ) +DD1 3 
DD22= (DD15-DD 1 7) / (30- (-30) ) * (DR- (-30) ) +DD 1 7 
DD3 1= (DD21-DD22) / (25~ (-25) ) * (DA- (-25) )+DD22 

c 

C CNO = DD3 1 IF MACH NUMBER =.9 

c 

CNO = D31 +(DD31-D31)ft(MACH-.6)/0.3 

c 

C cn=cnO+cnbftbetdeg+bo2vtft (cnr*r+cnpftp) 

c 

B 1 =BETDEG/20ftC 1 0XXB2+ (20-BETDEG) /20ftC 1 OXXBO 
B2=BETDEG/20ftC10NXB2+ (20-BETDEG) /20*C10NXB0 
B1 1= (B1-B2) / (25- (-25) ) * (DA- (-25) )+B2 
B3=BETDEG/20ftC10XNB2+ (20-BETDEG) /20ftC10XNB0 
B4=BETDEG/20*C10NNB2+ (20-BETDEG) /20ftC10NNB0 
B 1 2= (B3-B4) / (25- (-25) ) * (DA- (-25))+B4 
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d o o o o o r> 


B 1 3= (B 11 -B 1 2) / (30- (-30) ) * (DR- (-30) ) +B 1 2 


C 10=B 1 3 WHEN MACH NUMBER = .6 


IF (BETDEG.LT. -12.) BETDEG=- 1 2 . 

IF (BETDEG.GT. 12.) BETDEG= 12. 

BB1=BETDEG/12*C10XX2+(1 2-BETDEG) /12*C10XX0 
BB2=BETDEG/1 2*C 10NX2+(12-BETDEG)/1 2*C 10NX0 
BB1 1= (BB1-BB2) / (25- (-25) ) * (DA- (-25) )+BB2 
BB3=BETDEG/ 1 2*C 10XN2+ ( 1 2-BETDEG) / 1 2 * C 1 OXNO 
BBli=BETDEG/ 1 2*C 1 0NN2+ ( 1 2-BETDEG) / 1 2*C 1 ONNO 
BB 1 2= (BB3-BB4) / (25“ (-25 ) ) * (DA- (-25) ) +BB4 
BB 1 3= (BB 1 1 -BB 1 2) / (30- (-30) ) * (DR- (-30) ) +BB 1 2 


CIO = BB 13 IF MACH NUMBER = .9 


CIO = B 13 +(BB13-B13)*(MACH-.6)/0.3 

mMmmmmmmmmmttmmmmmmmmmm 

C end of interpolations 

c mmmmmmmmmmmmmmmmmmmm## 

RETURN 

END 


106 



APPENDIX D 

Computer Code for Simulation Comparison 
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== 

COMMOD_N 

THIS PROGRAM IS USED TO COMPUTE 

U 

. DOT, V. DOT, 

== 

== 


W. DOT, P 

. DOT, R. DOT USING THE HARV 

WIND-TUNNEL 

= = 

== 


MODEL. 




“ = 

bb 

INPUTS: 

FLIGHT DATA FILE 



=:= 

-= 


T 

TIME 



= = 

== 


MACH 

MACH NUMBER 



== 

== 


HAB 

ALTITUDE 



= = 

== 


QBAR 

DYNAMIC PRESSURE 



== 

== 


ALPDEG 

ANGEL OF ATTACK 




== 


BETADEG 

SIDESLIP ANGEL 



== 

== 


PH ID 

EULER BANK ANGEL 



== 

== 


THETAD 

EULER PITCH ANGEL 



== 

SB 


PS ID 

EULER YAW ANGEL 



= = 



P 

AIRCRAFT X-BODY AXIS ROLL 

RATE 

== 

== 


Q 

AIRCRAFT Y-BODY AXIS PI 

TCH RATE 

SB 

== 


R 

AIRCRAFT Z-BODY AXIS YAW 

RATE 

= = 

== 


DH 

STABILATOR DEFLECTION 



SB 

== 


DA 

AILERON DEFLECTION 



SB 

== 


DR 

RUDDER DEFLECTION 



BB 

== 

OUTPUTS: 





SB 

== 


UDT 

TIME DERIVATIVE OF 

U 


BB 

== 


VDT 

TIME DERIVATIVE OF 

V 


BB 

== 


WDT 

TIME DERIVATIVE OF 

W 


BB 

== 


PDT 

TIME DERIVATIVE OF 

P 


BB 

== 


QDT 

TIME DERIVATIVE OF 

Q 


= = 

== 


RDT 

TIME DERIVATIVE OF 

R 


BB 

II 

AUTHER 

Jl CHANG 

CAO 



:: 

== 


GRADUATE 

RESEARCH ASSISTANT 



BB 

== 


SCHOOL OF AEROSPACE ENGINEERING 



BB 

sss 


GEORGIA 

INSTITUTE OF TECHNOLOGY 



BB 

== 


ATLANTA. 

GEORGIA 30332 



= = 

== 

DATA 

JUNE, 1990 




== 


PROGRAM COMLAT 
PARAMETER (N=l 29) 

IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

REAL *8 MACH .MASS ,NZ ,IAS , HAB 

REAL *8 AU(N) ,AV(N) ,AW(N) ,AP(N) ,AQ(N) ,AR(N) 

REAL *8 APHID(N) , ATHETA (N) .APSID(N) , AALP (N) , ABET (N) 
REAL *8 TIME (N) , AMACH (N) , AHAB (N) , AQBAR (N) 

REAL *8 ADH(N) ,ADA(N) ,ADR(N) ,ADSB(N) 

REAL *8 UDT(N) ,VDT(N) ,WDT(N) ,PDT(N) ,QDT(N) ,RDT(N) 
real *8 atnetl (N) .ATNETR(N) ,lxe ,lye ,lze 
CHARACTER*80 HEADER 

CHARACTER AA 1 * 1 5 , AA2* 1 5 . AA3* 1 5 • A Ai** 1 5 . A A5* 1 5 . BB*50 
OPEN (UN I T=5 , F I LE= 1 MACH03 1 ) 

OPEN (UNIT= 6 ,FILE='DATA_NV) 


GENERAL CONSTANTS 


G =32. 17 1 * 

DTR =ACOS (- 1 .) / 1 80 . 


AIRCRAFT CONSTANTS 


CBAR =11.52 
MASS =1035.308 
XI = 23000 . 

Y I =151293. 
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Zl 

=169945. 

XZI 

=-2971. 

B 

=37.42 

CBAR 

= 1 1 .52 

S 

=400. 

XL 

=-3.56/12. 

YL 

=0. 

ZL 

=2.8/12. 

THZ 

=0.0 

LXE 

=-232.5/12. 

LYE 

=0. 

LZE 

=2-8/12. 

DT 

= 1 30. 

CSTAR 

=X 1 *Z 1 / (X 1 *Z 1 -XZ 1 StSt2) 

C41 ' 

=CSTAR*XZ 1 * (Z l+X 1 -Yl ) / (X 1 *Z 1 ) 

C42 

=CSTAR* (Z 1 * (Y 1 -Z 1 ) -XZ 1 **2) / (X 1 * Z 1 ) 

C43 

=CSTAR*XZI/XI 

C51 

= (Z 1 —X 1 ) /Y 1 

C52 

=XZI/YI 

C6l 

=CSTAR* (X 1 * (X 1 -Y 1 ) +XZ 1 **2) / (X 1 i-Z 1 ) 

C62 

=CSTAR*XZ 1 * (Y 1 -Z 1 -X 1 ) / (X 1 *Z 1 ) 

C63 

=CSTAR*XZI/ZI 

FLIGHT 

CONDITION 


READ (5,5) NUMBER 
DO 99 1=2,13 

READ (5,15) AA1 , AA2 ,AA3 ,AA4 ,AA5 
99 CONTINUE 

READ (5,25) BB 
5 FORMAT (14) 

15 FORMAT (5A15) 

25 FORMAT (A50) 


DO 400 1=1, N 

READ FLIGHT DATA FILE 

READ (5,10) T , MACH ,HAB , QBAR .ALPDEG .BETADEG ,PS!D .THETAD 
READ (5,10) PHID , PDEG , QDEG , RDEG ,VTOTAL,AX , AY ,GZ 
READ (5,10) DLADEG , DLADDEG , DLSDEG , DLHDD EG , DLNDEG , DLNDDEG , 

1 DLFDEG.DLFDDEG 

READ (5,20) DLRDEG , DLRDDEG .TNETL .TNETR 


: St St St St St St St St 

St Si Si Si s 

't St St Si Si St St Si Si Si St Si Si Si Si Si Si Si Si Si St Si St St Si St St St St Si Si St St 

T 

— 

TIME 

MACH 

— 

MACH NUMBER 

HAB 

— 

ALTITUDE 

QBAR 

— 

DYNAMIC PRESSURE 

ALPDEG 

— 

ANGEL OF ATTACK 

BETADEG 

— 

SIDESLIP ANGEL 

PS ID 

— 

EULER BANK ANGEL 

THETAO 

— 

EULER PITCH ANGEL 

PHID 

— 

EULER YAW ANGEL 

PDEG 

— 

AIRCRAFT X-BODY AXIS ROLL RATE 

QDEG 

— 

AIRCRAFT Y-BODY AXIS PITCH RATE 

RDEG 

— 

AIRCRAFT Z-BODY AXIS YAW RATE 

VTOTAL 

— 

VELOCITY 

DLADEG 

— 

AILERON DEFLECTION (AVERAGE) 

DLRDEG 

— 

RUDDER DEFLECTION (AVERAGE) 

DLSDEG 

— 

STABILATOR DEFLECTION (AVERAGE) 

TNETL 

— 

THRUST OF LEFT ENGINE 

TNETR 

— 

THRUST OF RIGHT ENGINE 

St St St St St Si St St St St Si St s 

: St St St Si St St St St St Si St St St St St St St St St St St St St St St St St St St St St St 


i'c St St St St St St St St St St St St St St 


'{ St St Si St St St Si Si Si St 


10 FORMAT (8E10.4) 
20 FORMAT (4E10.4) 
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Q_ O' £T 


u = VTOTALftCOS (ALPDEG*DTR) *COS (BETADEG*DTR) 

V = VTOTAL*S I N (BETADEG*DTR) 

W = VTOTAL*SIN (ALPDEGsVDTR) *COS (BETADEG^DTR) 


AU { 1 ) = U 
AV ( 1 ) = V 
AW ( 1 ) = W 

AP ( 1 ) = PDEG 

AQ ( 1 ) = QDEG 

AR ( 1 ) = RDEG 

APHID (1) 

=PH 1 D 

ATHETA(I) 

=THETAD 

APSID (1) 

=PS 1 D 

AALP (1) 

=ALPDEG 

ABET ( 1 ) 

=BETADEG 

TIME (1) 

=T 

AQBAR (I) 

=QBAR 

AMACH (1) 

=MACH 

AHAB (1) 

=HAB 

ADH (1) 

=DLSDEG 

ADA (1) 

=OLADEG 

ADR (1) 

=DLRDEG 

ADSB (1) 

=DLSBDEG 

ATNETL (1) 

=TNETL 

ATNETR(I) 

=TNETR 

CONTINUE 

CALL 1 AERO 

(HEADER) 

CALL 1 ENG 

(HEADER) 

DO 410 1=1 

. N 


ALPHA =AALP ( I ) 
BETA =ABET ( I ) 
MACH =AMACH ( I ) 
HAB =AHAB ( I ) 

DH =ADH ( I ) 

DA =ADA(l)/2. 
DR =ADR ( I ) 

DSB =ADSB ( I ) 
QBAR =AQBAR ( I ) 
TNETL =ATNETL ( I ) 
TNETR =ATNETR (I ) 


CALL F18M3 

1 

2 

3 

A 

5 


CDO 

, CYO , CLO 

,C10 

, CMO , CNO 

CLAD 

, CMAD , CLQ 

,CMQ 

*th ,tx 

CYB 

. CYR , CYP 

. C 1 B 

. C 1 R .Clp 

CNB 

. CNR , CNP 

, B02VT , C02VT. 

F ( 




ALPHA 

.BETA , MACH ,HAB 

f 

DH , 

DA , DR , 

DSB , 

DT ) 


U =AU ( I ) 

V =AV ( I ) 

W =AW ( I ) 

=AP (I) *DTR 
=AQ (I) *DTR 
=AR (I) *DTR 
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PH ID =APHID (I) *DTR 
THET =ATHETA (I) *DTR 
PS ID =APS I D ( I ) *OTR 

AA =ALPHA*DTR 

BT =BETA*DTR ■ '~V ' 

THX = (TNETL+TNETR) *COS (1 . 98*DTR) 

THY = TNETR*SIN (-1 ,98*DTR)+TNETL*SIN (1 .g8*DTR) 

FD = QBAR*S*C DO/MASS 

CB = QBAR*S*C02VT*CLAD/ (U*U+W*W) /MASS 

BQ = l+CB*COS (£A) *u-CB*COS (AA) *W* (CB*S I N (AA) *U) 

X / (1+CB*S I N (AA) *W) 

FU = R*V-Q*W-G*S I N (THET) ~FD*COS (AA) +THX/MASS 
X +QBAR*S*S I N (AA) * (CL0+C02VT*CLQ*Q) /MASS 

FW = Q*U-P*V+G*COS (THET) *COS (PH I D) -FD*S I N (AA) 

X -QBAR*S*COS (AA) * (CLO+C02VT*CLQ*Q) /MASS 

* WDT(I)- (FW+CB*COS(AA)*W*FU/(l+CB*SIN(AA)*W))+THZ/MASS 
UDT ( I ) = (FU+CB*S I N (AA) *U*WDT ( I ) ) / ( 1+CB*S I N (AA) *W) 

DALFA = (U*WDT ( I ) -W*UDT ( I ) ) / (U*U+W*W) 

CL = CL0+C02VT*(CLQ*Q+CLAD*DALFA) 

CM = CMO+C02VT*(CMQ*Q+CMAD*DALFA) 

Cl = C 1 0+C 1 B’’«BT+B02VT* (C 1 P*P+C 1 R*R) 

CY = CY0+CYB*BT+B02VT* (CYP*P+CYR*R) 

CN = CN0+CNB*BT+B02VT* (CNP*P+CNR*R) 

FL = QBAR*S*CL/MASS 

FX = -FD*COS (AA) +FL*S I N (AA) 

FY = QBAR*S*CY/MASS 

FZ = -FD*SIN(AA)-FL*COS(AA) 

FP = QBAR*S*B*C1/XI+MASS*(YL*FZ-ZL*FY)/XI 

FQ = QBAR*S*CBAR*CM/YI+MASS* (ZL*FX-XL*FZ) /Yl 

FR = QBAR*S*B*CN/ZI+MASS*(XL*FY-YL*FX) /z I 

VDT(I) = P*W-R*U+G*COS (THET) *S IN (PH I D) +F Y+THY/MASS 

YY1 = CAl*P*Q+C42*Q*R+Ci*3*FR+CSTAR*FP+Cit3>VLXE*THY/ZI 
1 -Ci»3*LYE*THX/ZI 

1 +CSTAR* (LYE*THZ-LZE*THY) /X I 

YY2 = C5 1 *P*R+C52* (R*R-P*P) +FQ+ (LZE*THX-LXE*THZ) /Y I 
YY3 = C6 1 *P*Q+C62*Q*R+c63*FP+CSTAR*FR+C63/X I * (LYE*THZ-LZE*THY) 
1 +CSTAR/ZI *LXE*THY 

1 -cstar*lye*thx/zi 

PDT(I) =YY1/DTR 
QDT(I) =YY2/DTR 
RDT(I) =YY3/DTR 

b 10 CONTINUE 

M=N-2 

DO 710 1=1, M 

WRITE (6,910) TIME (I) , UDT ( I ) , VDT ( I ) ,WDT(I) ,PDT(I) , QDT ( I ) , 

1 RDT(I) 

710 CONTINUE 
910 FORMAT (7 E 10.4) 

STOP 

END 
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APPENDIX E 


Comparison of Cmo(t) for Run 5, Mach 9 Flight Trajectory 


The piloted simulated maneuvers comparison of Chapter 10 shows that the angular pitch 
accelerations from the wind-tunnel data and the analytical model have about the same shape 
but at times have a fairly large distance between them; this is due to a small error in fit being 
multiplied by a large dynamic pressure at Mach 0.9. We found that the largest differences 
are equivalent to approximately a half degree change in stabilator deflection. In this 
appendix we present some details from Run 5 which is a turn reversal maneuver to show 
that the differences are due to a small difference of about 0.006 or less in the values of 
Cm 0 (t)- As can be seen from the modeling fits shown in Figures 5. 1-5.6 modeling errors 

of this magnitude are present in C mo at all Mach numbers. 

The time history of C m is p otted in Figure E.l(a) for Run 5 showing maximum 
differences of about 0.01 magnitude. Removing the effect of dynamic pressure we note 
that the maximum difference in is about 0.006 as shown in Figure E.l(b). The time 
histories of the angle of attack and the stabilator angle values are presented in Figures E.l 
(c) and (d), respectively. The angle of attack has values between 1 and 5 degrees and the 
stabilator angle has values between 1 degree and -3.0 degrees. Consequently, the 
analytical models of Chapter 5 governing Run 5 are CM0X29(a), CM0X9(a) and 
CM0N59(a) presented in Figures 5. 5,5.6 and 5.8. 

As can be seen from Figures E.l(b) and E.l(c) the small difference 0.006 can be 
made up by small changes in the stabilator angle. Such a small difference in the stabilator 
angle would have negligible bearing on analysis study results using the analytical models as 
compared to those obtained using the wind-tunnel data. 
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ALPHA 


Comparison of C rn() (t): Turn Reversal Maneuver @ M=0.9 (Run 5 , 10 October 1987) 



Wind-Tunnel Analytical 

Figure E.l: Comparison of C^t): Turn Reversal Maneuver @ M=0.9 (Run 5, 10 October 1987) 
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